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ABSTRACT 


We describe a model of pulsar magnetospheres which has evolved induc- 
tively from the work of Sturrock, where the radiation is produced near 
the surface of a neutron star. Those observations which are central to 
our arguments are briefly reviewed. As a foundation for our work we 
discuss some of the theoretical ideas of others, particularly those of 
Sturrock . 


The braking index n and period-pulse-width distribution of pulsars 

are first re-investigated by relaxing the conventional assumption that 

Ry = R^, where Ry is the radius of the neutral points marking the 

transition from closed to open magnetic field lines, and R^ is the 

radius of the light cylinder. This is replaced by the parametrization 
1 — Tl 

Ry = R^ Rl , where R * is the neutron star radius. If the radio fre- 
quency radiation is created near the surface and beamed along open field 
lines, it is found that a good fit to the period-pulse-width distribution 
can be obtained for 71 in the range 0.5 5 7] ^ 0.7- The relation n = 
1+27] then gives n = 2.2 ± 0.2, which is in good agreement with the 
values measured for the Crab pulsar. 

Recent observational studies concerning timing irregularities and 
the dispersion measure of the Crab pulsar and wisp motions in the Crab 
nebula indicate that pulsar magnetospheres may contain large amounts of 
non-relativistic gas. We point out that substantial accumulation can 
occur only in the vicinity of the force-balance radius R^, where 
gravitational and centrifugal forces cancel. It is argued that Ry » 

For this model, the braking index n = 7/3 and the expected period- 
pulse-width distribution and period-age distribution agree well with 
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observational data. A glitch is interpreted as an instability leading 
to the ejection and/or dumping of accumulated gas. 

We have constructed approximate models for the magnetic field struc- 
ture, considering in detail the cases of aligned and orthogonal spin and 
magnetic axes. If the magnetic field is initially closed (out to the 
light cylinder) gas will collect where there is "microscopic stability" . 
Such locations exist for r > R^. Most of this region is "macroscopically 
unstable" since, when sufficient gas is collected, the centrifugal force 
overcomes the magnetic stress and converts closed field lines into open 
field lines. However, there is a small region near r = R^, both 
microscopically and macroscopically stable, where gas can collect and 
remain . 

Such a magneto spheric structure leads to a braking index n = 7/3* 
in good agreement with the Crab pulsar. We calculate the polar cap 
boundaries, and the resulting period-pulse-width distribution agrees well 
with an improved sample of pulsars. The surface magnetic field strengths 
deduced from this model are roughly a factor of ten less than pre- 
vious estimates; for the Crab pulsar we find ey 2 X 10 11 gauss. 

In conclusion we discuss the ability of the resulting model to 
account for the observed features of pulsars in general and the Crab 
pulsar in particular. We mention some shortcomings, and suggest possible 
future avenues of investigation. 
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Chapter 1 


INTRODUCTION 

Ever since Gold (1968) made the cogent suggestion that pulsars are 
basically rotating neutron stars, theoretical studies have followed two 
more or less distinct lines. First are those which attack the problem 
of a rotating, conducting, magnetized sphere from first principles. The 
aim of these papers has been to discover the structure of the electro- 
magnetic fields and plasma surrounding such a neutron star (the magneto- 
sphere) . It is assumed that once this has been done, the mechanisms 
which produce the observed radiation will be discemable. 

The second class of studies contains those which examine specific 
processes which might produce radio, optical, x-ray, or gamma radiation 
in the vicinity of a neutron star. The magnetic field structure is often 
prescribed by fiat , and may or may not be related to model magnetospheres 
studied by others. Broadly speaking, radiation mechanisms may be divided 
into two groups, depending on whether the radiation is thought to arise 
near the surface of the star or at the light cylinder. 

In this thesis we describe a model of pulsar magnetospheres which 
has been developed inductively around the hypothesis that radiation is 
produced near the surface of the star, specifically, at the magnetic poles 
(Radhakrishnan et al . 1969 ; Komesaroff 1970; Sturrock 1970, 1971a) . 
Analysis of various data has led us to modify our model in stages, to 
the point where it can account for a number of important observations in 
a unified manner. 

A short review of the relevant observations is given in Chapter 2. 

To provide a background and in order to introduce some of the ideas and 
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definitions which will be used frequently, we discuss briefly some of 
the important theoretical work of others. The model of Sturrock (I97O, 
1971a, 1971b) is described in some detail, since it provides the frame- 
work upon which this work is based. 

The bulk of this thesis consists of three papers by the author and 
his advisor, Professor Peter Sturrock. They are presented essentially 
as published, with several appendices to elaborate upon certain calcu- 
lational details . 

In Chapter the braking index n and the period-pulse-width (W-T) 
distribution of pulsars are re-investigated by relaxing the conventional 
assumption that R^ = R^, where R^ is the radius of the Y-type neutral 
points marking the transition from closed to open magnetic field lines, 
and R^ is the radius of the light cylinder. This assumption is re- 

T| 1-n 

placed by the parametrization R^ = R^ R # where R^ is the neutron 
star radius. Comparison is made with the observed braking index of the 
Crab pulsar, the W-T distributions of two samples of pulsars, and with 
the number of pulsars showing an interpulse. It is found that a value 
of T1 in the range 0.5 < 71 ^ 0.7 is consistent with all of the data. 
(This material is from Paper I, Roberts and Sturrock 1972 a.) 

The magnetospheric structure discussed in Chapter 3 (where R^ « V 
is radically different from that usually assumed (R^ = R^) . In Chapter 4 
it is shown that such a configuration is a natural consequence of the 
presence in pulsar magnetospheres of the large amounts of non-relativistic 
gas which were inferred from observations by Scargle and Pacini (I971). 

It is argued that in such a case Ry sa R^g, where is the force- 

balance distance, the radius at which the gravitational and centrifugal 
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forces on a corotating particle cancel out. For this model n = 7/3 j 
and the resulting W-T and period-age distributions agree well with 
observational data. A glitch is interpreted as an instability leading 
to the ejection and/or dumping of accumulated gas. (Paper II, Roberts 
and Sturrock 1972b . ) 

This model is explored in greater detail in Chapter 5 . We have 
constructed approximate models for the magnetic field structure, con- 
sidering in particular the aligned and orthogonal cases. The stability of 
single particles ("microscopic stability") is studied in detail. The 
more difficult question of the stability of large amounts of collected 
material ("macroscopic stability") is discussed in simple physical terms. 
We find the consequent magnetic field structure and braking index 
(n = 7/3) and the boundaries of the polar caps. The resulting W-T 
distribution and number of pulsars expected to show interpulses are com- 
pared with an improved sample of pulsars. (Paper III, Roberts and 
Sturrock 1973*) 

In Chapter 6 we review our modifications of Sturrock' s model, 
commenting on the limitations as well as the strengths of the results. 

We discuss areas in which our work may be extended . Finally, we comment 
briefly on those observations, as yet unexplained, which we believe to be 
of crucial importance to a fuller understanding of the physics of pulsars. 
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Chapter 2 


REVIEW OF OBSERVATIONS AND PREVIOUS THEORETICAL WORK 
I . Observations 

There is certainly no dearth of observational facts to challenge the 
pulsar theorist. Below we present a brief review of those which are 
especially relevant to our work. The wealth of data which are not pre- 
sented have been more than amply reviewed by Hewish (1970), Huguenin, 
Manchester, and Taylor (1971), Ginzburg (1971), ter Haar (1972), and 
Backer (1973), among others. 

Each of the 89 pulsars presently known was first detected by its 
extremely periodic bursts of radio frequency radiation. The observed 
repetition period T of the pulses lies between 0.033 sec (Crab, PSR 
0531+21) and 3-75 sec (PSR 0525+21), with a median of 0.655 sec. In all 
cases where the time derivative T has been measured it is positive, and 
the "age" t = T/T varies between 2.5 X lO " 5 yr (Crab) and 2.5 X 10 yr 
(PSR 0809t7l+) . 

Another observational fact is the relatively short duration W of 
the pulses compared to the period (W/T » 0.05). Six of the pulsars (in- 
cluding the Crab pulsar) exhibit an "interpulse" — a second, smaller 
pulse which comes roughly half-way between consecutive main pulses. It 
is also possible to classify pulsars according to pulse shape, and whether 
they show the "drifting sub-pulse" phenomenon, in which features of the 
pulse move through the overall envelope from pulse to pulse (Taylor and 
Huguenin 1971 ; Backer 1973) . 

A feature of the radio emission unique to the Crab pulsar is the 
"precursor". It occurs about 1.6 msec before the main pulse, has no 
optical counterpart, and is observed only below 600 MHz. Below 400 MHz 
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its energy is comparable to that of the main pulse, and its duration is 


roughly five times longer. The polarization is almost 100 percent linear, 
and the position angle is the same as that of the main pulse. 

The extremely accurate timing data which has been accumulated for 
the Crab and Vela (PSR 0833 - ^5) pulsars exhibit additional features which 
require explanation. First, the "braking index" n (n = axiVcn = 2-TT/T ; , 
which measures the frequency dependence of the torque which slows the 
neutrons star's rotation (torque ce a/ 1 ) , has been measured for the Crab 
pulsar; it is found to have a value n = 2.4 ± 0.2 (Boynton et al . 1972 ) . 
Second, the slow-down is not smooth, but exhibits both seemingly random 
"noise" (Crab) and larger, apparently discontinuous, jumps in frequency. 
Two such "glitches" have been observed for each of these pulsars; 

-Q -6 

Aou/cd ss +10 for the Crab pulsar, and Aoc/cn » +10 for the Vela pulsar. 
According to Boynton et al . (1972), the noise seen in the Crab may be 
well represented by mini-glitches of AuV'cn s 2 X 10 ^ occurring at least 
once per week . 

Differences in the times of arrival of a pulse at two different 
frequencies may be expressed in terms of the dispersion measure 


t(v i )-t(v j ) 



( 2 . 1 ) 


When the radio frequency is much greater than the plasma frequency 

f 2/ vl/2 
v = (n e / jtm ) , 

p v e e ' ’ 


D. . 

1J 



10~2 . 88 ^ d £ _ 


( 2 . 2 ) 
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Careful observations of the Crab pulsar (Rankin, Counselman, and Richards 

1971 5 Rankin, Campbell, and Counselman I 97 I ; Rankin and Counselman 1973) 

1 7 -p 

show that its dispersion measure (nominally about 2.36 x 10 sec Hz , 

so that In dJt sa 56.8 parsec cm is not strictly constant, but shows 

k 

variations on the order of 2 parts in 10 on time scales of months . One 
such change in integrated plasma density occurred during a two month 
period beginning around the time of the September 1969 glitch. However, 
similar fluctuations have occurred without being associated with timing 
irregularities, and no increase in dispersion measure followed the Crab 
glitch of October 1971 ■ 

Careful observations of small-scale features ("wisps") in the Crab 

nebula near the pulsar have revealed motions and brightenings that may be 

associated with timing irregularities of the pulsar (Scargle 1969 ; Scargle 

and Harlan 1970; Scargle and Pacini 1971) . Correlations are claimed for 

small irregularities occurring near March 10 and June 1, 1970, as well as 

with the large glitch of September 1969 * but there was no activity following 

the large glitch of October I 971 . The energy required for these events 

4l 

is claimed to be at least 10 ergs (Scargle and Harlan 1970) . 

The high brightness temperature of the rf radiation was seen from the 
outset to require a coherent mechanism. (Assuming the radiation comes 
from a region of linear dimension less than cW, the brightness temperature 
is about 10 2 °-10 2 ^ °K . ) It appears that the rf signals display (in 
addition to interstellar scintillation effects) fluctuations intrinsic 
to the source, notable examples of which are the drifting subpulses, the 
"nulling" phenomenon (Backer 1970 ) , and the "giant pulses" of the Crab 
pulsar (Sutton, Staelin, and Price 1971 ) . There is a great deal of 
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spectral information concerning pulsars. In general, there is a low 
frequency cut-off around 100 MHz and the radio flux is a decreasing 
function of frequency. In addition, the rf radiation is frequently 
strongly linearly polarized, and is often found that the electric vector 
sweeps through a large angle during a pulse ( Radhakrishnan et al . 1969) . 

Pulsed radiation in other frequency bands, coincident (after correc- 
tion for dispersion delays) with the rf pulses, has been detected from 
the two fastest pulsars (Crab and Vela) . While the Vela pulsar has been 
detected only in x-rays (Harnden et al . I972), the Crab has been observed 
in the optical (Cocke et al . 1969), x-ray (Fritz et al. 1969) and gamma- 
ray (Hillier et al . 1970 ) spectral ranges. The ratio of energy in the 
interpulse to that in the main pulse, which is small in the radio band, 
grows steadily from the optical through the gamma-ray, where it is larger 
than unity. 

The optical pulse of the Crab pulsar is exceedingly stable. The main 
pulse-interpulse configuration is asymmetric, and exhibits a cusp-like 
main peak, unresolved on a time scale of 30 microseconds (Papaliolios 
et al . 1970 ). The optical pulse shows some linear but no circular polari- 
zation. The electric vector sweeps through an angle of 130 ° during the 
main pulse, and is parallel to that of the rf pulse at the peak (Manchester 
1971a) . The spectrum of the optical pulses has a broad maximum in the 
visual part of the spectrum, and turns over in the infrared. The spectrum 
apparently extends into the x-ray band and from there to the gamma-ray 
band. There was no change in the optical pulse shape or intensity during 
or following the September 1969 glitch. 
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[I. Pulsar Electrodynamics-Canonical Models 

One possible approach was taken by Pacini (1968) and Ostriker and 
Gunn (1969; Gunn and Ostriker 1969, 1970 ). They examined the dynamics 
of an object possessing gravitational or magnetic fields not symmetric 
about the rotation axis. In a vacuum, electromagnetic and gravitational 
energy is radiated at rates given by 

dE md /dt = " § “j 2 /c3 ( 2 - 3 ) 

and 

dE gq /dt = " G D 2 ^/ c5 • (2 .h) 

Here m and D are the components of the magnetic-dipole and mass- 
x x 

quadrapole moments perpendicular to the rotation axis. Such an object, if 
rotating almost rigidly, will also loose angular momentum at a rate given 
by 

dJ/dt = <d 1 dE/dt . ( 2 . 5 ) 

Thus the braking indices for these models are n = 3 (magnetic dipole 
radiation) and n = 5 (gravitational quadrapole radiation). 

Surface magnetic fields of 100 gauss are not unreasonable for young, 
massive stars, the likely progenetors of neutron stars. Flux conservation 

in a collapse from a radius of 10^ cm would produce a surface field 

12 6 

B„ on the order of 10 gauss at R„ =10 cm. Thus the magnetic moment 
R 3 b„ fej 10 30 emu for the resulting neutron star, and Eq. ( 2 . 3 ) becomes 

dE^/dt bs -10 2 ®' 3 co* erg s 1 . (2.6) 
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p o -Crab *37 *5 -1 

The Crab pulsar has a) 10 ‘ J s , so E , & -10 J erg s . This 

ma 

interpretation is strengthened by the observation (Gold 1969 ) that if it 
is slowing down, it is losing rotational energy at the rate 


d 2 • -1 

- -j-jr [(1/2) I^cd ] = -I^cd cd erg s 


(2.7) 


kb ’i 2 

Neutron star models have moments of inertia I =- 10 g ci , the Crab 

. -Q 6 -2 37.2 

is observed to have cd -10 s , so the energy loss is « 10 

i 

erg s , . In addition, this is about the right amount of energy that is 
required to power the Crab nebula (Shklovsky 1968; Haymes et al . 1968; 
Wang and Sartori 1973) • 

Combining Eqs . (2.3) and (2.7)> we may express the surface magnetic 
field in terms of the observed slow-down rate: 



It is rather remarkable that although T varies over three orders of 

magnitude and T over two orders, the magnetic field (evaluated for a 

given neutron star model from Eq . [2.8]) varies by only a factor of four 

1/2 3 

for the 25 measured pulsars. However, the large range of /R , 

especially for small neutron stars masses, makes it impossible to deter- 
mine magnetic fields to better than 3 -^- orders of magnitude (Greenstein 
1972) . 

Gravitational radiation will always dominate over magnetic radiation 
if the frequency is sufficiently high. However, the measured braking 
index of n w 2,4 for the Crab pulsar implies that this is not important 
at the present epoch, and we will not discuss it further. 
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The simplest model of a pulsar magnetosphere starts with a neutron 


star possessing a dipolar magnetic field, its magnetic axis parallel 
(or antiparallel) to its spin axis. Although it is obvious that such 
a time-independent axisymmetric model cannot account for the pulses, it 
is usually assumed (and Mestel [I97I] has shown this to be very plausible) 
than an oblique rotator will be qualitatively the same as an aligned 
rotator . 

The pioneering calculation was that of Goldreich and Julian (1969)* 
who noted that the extremely good conductivity inside the star requires 
that 



If there are no charges surrounding the star, the potential for r > 
may be found from Laplace's equation. The resulting electric field is 
quadrapolar, there is a surface charge 

B*R*03 

CT = - "i t - ■ cos 0 , (2.10) 


and outside the star 


03 R 8 B 2 

-i —* * X 

E • B = = — cos 0 . 

c r 


( 2 . 11 ) 
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However, for the magnetic fields ( B^ 10 gauss) in pulsars such a con- 
figuration is highly unstable, since the non-zero electric force along 
the magnetic field is enormously greater than the gravitational force. 
Unless quantum effects at the surface are sufficient to bind the charges 
(Ruderman 1971 ; Ginzburg and Usov 1972), they are pulled off into the 
magnetosphere and behave as a conducting extension of the star. In this 
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case, as long as inertial and gravitational forces are small compared with 
electric forces, the Lorentz force law and infinite conductivity require 

—4 “4 

that E ■ B = 0 outside the star as well as inside. The resulting motion 
of charged particles is a combination of corotation and motion along field 
lines. When the contribution of the corotating charge to the magnetic 
field is included, the charge density outside the star is approximately 


n -n 


+ 


-Tn • B 
2itec 


. l-( cur/c )^sin^8 J 


( 2 . 12 ) 


This applies only to the corotating part of the magnetosphere, which con- 
sists of those magnetic field lines which lie entirely within the "light 
cylinder" . The radius of the light cylinder, where the co rotation speed 
is equal to the speed of light, is given by 


= c/cd = cT/2jt . (2.13) 

Field lines which penetrate the light cylinder must be open in order that 
particle velocities remain less than c. 

The torque N on the star is obtained by integrating the Maxwell 
stress tensor over a spherical surface centered on the star. The result 
is 

X “ b # \ 6 0) 3 /8c 3 = £ [n,( W 3 ] 2 V* (2 ' l4) 

where is the radius of the neutron star, and the polar magnetic 

field strength at the surface. From Eq. (2.14) we see that the braking 
index in the Goldreich-Julian model is n = 3 . Equation (2.14) leads to 
an energy loss comparable to that from a vacuum rotator. 
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Although this model surely possesses some of the most important 


features, it is not a self-consistent solution of Maxwell's equations 
and the dynamical laws. Attempts to extend this analysis have been made 
by Scharlemann and Wagoner (l973)j Michel ( 1973 a, 1 973 b ) and Julian (1973), 
who also neglect inertial and gravitational forces. However, these calcu- 
lations seem to require specification of the current at the surface of the 
star or the configuration of the magnetic field at infinity. Thus it is 
not clear that a unique solution exists, at least as the problem has been 
posed. (it is contended that inclusion of inertia will make the solution 
unique, but this has yet to be demonstrated.) 

The equation of motion 

d(l^o))/dt = -N = -An / 1 (2.15) 

may be integrated easily. Assuming that A and 1^ are independent of 
time, and the present epoch to be t = 0 , 

m(t) = cn(0) [1 + (n-l)t/T] 1// ^ 1 , ( 2 . 1 6 ) 

where 

T = T/fj = — CD/ CD j . (2.17) 

t =0 t =0 

For magnetic braking (n = 3)> the frequency becomes singular at a time 
t = -t/ 2. For the Crab pulsar, this agrees reasonably well with the 
birth of the Crab nebula in 1054 AD. Or taking the age of the Crab pulsar 
as 9^-9 years and n = 3> co(birth) = 2.3 co(0). Lowering n towards 
the observed value decreases oo( birth), but not very much. However, the 
largest n that could have operated from birth is about 3*4. 
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An excellent review of pulsar electrodynamics and pulsar dynamics 
has been provided by Ruderman (I972). 

Ill . Radiation Mechanisms 

Although these two models explain reasonably well the kinematics of 
pulsar slow-down and the energy budget of the Crab nebula, neither 
attempts to account for the observed pulses. Many theoretical ideas con- 
cerning the nature of the radiation mechanism were advanced soon after 
the discovery of pulsars (see, for instance, Hewish I97O) . Most of these 
ideas have died quietly, either from neglect or under the weight of ad- 
verse evidence. One school of thought, stimulated by an early suggestion 
by Gold (1968), holds that radiation is produced by relativistic particles 
in the vicinity of the light cylinder. This difficult problem has been 
attacked by Shklovsky (I970), Eastlund (1968, I97O), Lerche (1970a, 1970b), 
Pacini and Rees (1970), and others. However, Michel (1973a) has recently 
shown that in a self-consistent inertia-less model where all the plasma 

corotates and does not flow along field lines, the radiation arising at the 

8 

light cylinder is too low by a factor of about 10 to be the source of 
pulsar rf radiation. 

Proceeding in the direction first indicated by Radhakrishnan et al . 
(1969)5 other theorists (Komesaroff 1970; Pacini and Rees I97O; Sturrock 
1970, 19Tla; Tademaru 1971) have explored in detail the idea that the 
short duty cycle and linear polarization sweep of the rf radiation are 
consequences of emission near the neutron star surface at the magnetic 
poles. These models require that the flow of charged particles at the 
polar caps should be bunched. It has also been proposed (Chiu and Canuto 
1971 ; Ginzburg et al . 1969) that an alternative way of obtaining radiation 
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of high brightness temperature is through some kind of maser mechanism. 

ter Haar (1972) has reviewed these and other suggestions, and the 
reader is referred to his article for a thorough discussion of the merits 
and shortcomings of these and other proposals. In what follows we shall 
be concerned with those models in which the radiation is produced close 
to the surface of the neutron star at the magnetic polar caps. 


IV. Sturrock's Pulsar Model 

Sturrock (I 97 O, 1971b) has constructed a model of pulsars 

extending the work of Goldreich and Julian ( 1971)7 but relaxing their 
assumption that the electric and magnetic fields are everywhere orthogonal, 
In the Goldreich- Julian model there is a toroidal (B ) component so that 


f 5 


the integral $ B • ds around the light cylinder is non-zero. This means 
that curl B is non-zero, which implies that there are currents flowing- 
in the system. The distribution of these currents, which flow from the 
magnetic polar caps out to infinity, may be roughly inferred from the 
magnetic field configuration. There must be a positive current from an 
angular ring and a negative current from the center of the polar cap, or 
vice versa, depending on whether the magnetic and spin axes are aligned 
or anti-aligned. This view was advanced also by Goldreich and Julian. 

At this point Sturrock departs from Goldreich and Julian by noting 

that a laboratory device such as an electron gun will not give a current 

unless a voltage is applied to it . In fact there is a source of voltage 

at the polar cap, namely the rotation-induced electric field. For the 

case of the Crab pulsar, ions would accelerate to about 10 eV per 

electronic charge. Because of radiation reaction, electrons will reach 
lip 

only about 10 eV. This acceleration occurs very close to the star, 



probably within a kilometer of the surface. Since the ions and electrons 


move along the curved magnetic field lines, they are copious emitters of 

12 

high energy gamma-rays (of order 10 eV) . 

The key mechanism in Sturrock's model is the conversion of these high 

energy gamma rays into electron-positron pairs in the intense magnetic 

field. The gamma rays are produced by curvature radiation and initially 

move parallel to the magnetic field. However, since the field is curved, 

they soon see a transverse component of magnetic field. When the condi- 

tion E (eV) B (gauss) £ 10 * is met, the gamma ray will annihilate 

y ± 

into an electron-positron pair. 

There are two important consequences of this pair creation. One is 

that the electrons and positrons are produced with non-zero pitch angles, 

and will therefore radiate by the synchrotron mechanism, possibly in the 

gamma, x-ray, or optical regions of the spectrum. Some of the gamma rays 

■f “* 

produced by the secondary electrons will in turn annihilate into e e 

7 5 

pairs, and so on. The resulting cascade produces about 10 y secondaries 
per primary electron, and the total number of relativistic leptons in- 

41 -1 

jected into the Crab nebula is around 10 s . This agrees with the number 
required to produce the synchrotron radiation from the nebula (Shklovsky 

1968) . 

The other consequence follows from the fact that some positrons will 
probably be turned around by the electric field and flow toward the sur- 
face. If this happens there will be a two-stream configuration which is 
known to be very unstable. This is likely to lead to bunched rather than 
steady flow. If bunches of electrons are moving along curved field lines, 
then radio emission will be produced by the mechanism outlined by 
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Radhakrishnan and Cooke (1969), and more fully analyzed by Tademaru (1971). 
This radiation is beamed along the direction of the magnetic field, and 
is polarized along the projection of the magnetic field seen by the obser- 
ver. Thus, as pointed out by Radhakrishnan et al . (1969), the electric 
vector will sweep in position angle as the line of sight of the observer 
passes through the polar cap. 

The hypothesis that pair creation is responsible for the radio 
emission may be tested (Sturrock 1971b) . Expressing the magnetic field 
and the electron energies (via the electric field) in terms of the period 
T and the age r, the condition that pairs can be created in the elec- 
tron polar zone (EPZ) may be written as 


* < *f 3 ^ 3 


(2.18) 


Similarly, for pair creation in the ion polar zone (IPZ), the requirement 
is 


< T iPC 


= 10 


-17.3 t 2/3 -2/3 




(2.19) 


2 

Here 1^ (g cm ) is the moment of inertia of the neutron star. These 
conditions are compared with the observed period-age distribution in 
Figure 2.1. It is seen that radio emission ceases approximately where one 
would expect pair production to end in the EPZ. The condition (2.19) is 
satisfied for only the Crab pulsar, and Sturrock (1970) has suggested that 
this offers an explanation of why the Crab pulsar alone shows a precursor 
(a similar suggestion has been advanced by Pacini [I97I]). 

The spectrum computed on the basis of this model account fairly well 
for the radio emission. The x-ray emission of the Crab pulsar may be 
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understood reasonably well if it comes from the IPZ . 


On the basis of this model, Sturrock made three predictions. First, 
that gamma radiation should be detectable from the Crab pulsar. This has 
been verified. Second, that x-rays should be observable from the Vela 
pulsar. This, too, has been verified. The third prediction was that the 
polarization of the optical radiation should have an electric vector 
orthogonal to that of the radio emission. This has been found to be 
untrue (Manchester 1971 a ) > and has led to the abandonment of the mecha- 
nism originally proposed for the optical radiation. Recently, Turk, 
Sturrock, and Petrosian (1973) have shown that coherent curvature radia- 
tion from first or second generation secondary electrons can account for 
the optical radiation from the Crab pulsar. 

However, Sturrock' s model predicts a braking index of n = 3 , and 
pulse widths that are much smaller than those which are observed. The 
way in which we have modified his model to account for these observations, 
as well as provide a basis for understanding the irregularities in the 
timing data, is described in the subsequent chapters. 
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FIGURE CAPTION 


Figure 2.1. Distribution of pulsars in period T and age r, compared 
with the predictions of Eqs. (2.1) and (2.2) (Sturrock 
1971b). ^,S-type structure; /S., C-type structure;Q ; D-type 
structure, in the notation of Taylor and Huguenin (1971) . 
The theoretical curves are based on the assumption that 

R Y= V 
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Figure 2.1 



Chapter 3 


THE BRAKING INDEX AND PERIOD-PULSE-WIDTH DISTRIBUTION 

I . Introduction 

Although there is some uncertainty concerning the braking index for 

the Crab Pulsar, it is generally agreed that the observations cannot be 

reconciled with any existing theory (Richards et al . 1970; Duthie and 

Murdin 1971 > Roberts and Richards 1971) • Likewise, there is a problem 

concerning the distribution of pulsars in period (T, sec.) and pulse- 

width (W, sec.). As noted by Gunn and Ostriker (1970), one expects that 

1/2 

the observations should satisfy the law W a T if the (RF) radiation 
is produced near the polar caps defined as those areas of the surface of 
a neutron star coupled by magnetic-field lines to the light cylinder. 1 
The purpose of this chapter is to show that both sets of observations 
may be understood in terms of a modified polar cap model. 

II . Model — Braking Index 

It is generally assumed that the transition from closed to open 
magnetic field lines occurs at the light cylinder, as proposed by Goldreich 
and Julian (1969). H is clear that the radius R^. of the "Y-type" 
neutral points (Sturrock and Smith 1968) , which characterize the transi- 
tion from closed to open field lines, cannot occur outside the light 
cylinder, which has radius R^ = cT/2rt. However, it is not at all clear 
that Ry cannot be less than R^ . Sturrock (1971b) has suggested that 
one consider to be determined as follows: 


1 0 7 

Gunn and Ostriker find that a better fit is W « T ‘ , which corres- 
ponds to T| = 0.6 in the model to be described. 
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(3.1) 


R - R 1_T ) R ^ 

Ry = R * ^ , 


where T\ is a parameter in the range 0 < T| ^ 1 . If one notes that the 

-3 

magnetic field strength varies as r in the range R < r < R and 

7T Y 

~p 

as r in the range Ry < r < R^, and if one notes that one should 
retain the assumption that the magnetic-field spiral-angle is about 45° 
at the light cylinder, one finds that the formula for the magnetic 
torque N given by Sturrock (1971a) is replaced by 


n - 1 B l 4 v 2 1’ 1 


(3.2) 


where R^(cm) is the radius of the star and B^( gauss) is the value of 
the magnetic field at the surface of the star. On using equation (3.l)> 
this becomes 


12 iw-21) R -!-2T1 
N = 2 B * R * \ ' 


(3.3) 


-i-2n 

which shows that N a T . Hence the braking index n is given by 


n = 1 + 2T) 


(3.4) 


III . Period-Pulse-Width Distribution 

We next consider the period-pulse-width distribution to be expected 
on the basis of equation (3.1), by extending the analysis of Henry and 
Paik (1969). We assume that radiation is due to charge sheets moving 
along open magnetic-field lines at the polar caps. We consider, for 
definiteness, that the dominant contribution is produced close to the 
surface of the star, but the conclusions are not greatly affected if the 
radiation is produced within a few star- radii. 
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If the magnetic dipole axis makes an angle © with the rotation 

axis, field lines which are tangent to a cylinder of radius R^, parallel 

to the rotation axis, define a polar cap at each magnetic pole. The semi- 

1/2 - 1/2 

major-axis is of angular size R^ R^ and is in the direction of in- 
creasing longitude (referred to the rotational axis); the semi-minor- 

1/2 - 1/2 

axis is of angular size R^ G( ®) and is in the direction of 

increasing latitude. G(0) varies slowly over the range G(0) =1 to 
G(jt/2) = .620. If the radiation polar diagram is defined by the cone 
which is tangent to the field lines passing through the boundaries of each 
polar cap, the major and minor half-angles of this cone are 


3 „l/2 

R * 


te, = i r: ' " R 


1/2 

■Y ' 


♦b = I R i /2 V 1/2 ° ( ®) 


(3.5) 


The details of this calculation are given in Appendix A. 

If <j) is the angle between the rotation axis and the line of sight, 
the pulse-width W is given by 


W 


T te, 

rt sin () 


1 


_ (<t ~ ©V 


ril/2 



(3.6) 


for one polar cap, and a similar expression (©-*«-©) for the other 
polar cap. (if W is found to be imaginary, radiation is not detectable 
by the observer.) 

If we assume that the rotation axis is randomly oriented, and if we 
assume that the magnetic axis is randomly oriented with respect to the 
rotation axis, we find (neglecting certain terms of order W/T) that 
the probability that radiation from either pole is observed and has a 
pulse width in the range W to W + dW is 

P w (W,T)dW = T -1 L(x) dW . (3.7) 
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In this expression. 


x = ky 


where 



and 



The function L(x) is given by 


L(x) 



d® 


sm 


0 G( ®) 


, . 2_v 1/2 

( x-sm ®) 


) 


where 


1/2 

® = arc sin x , x ^ 1 

m ’ ^ 

= tc/2 , x > 1 . 


(3-8) 


(3-9) 

( 3 . 10 ) 


(3.ii) 


( 3 - 12 ) 


In order to make comparison with observational data, it is easier 
to consider the probability that a pulsar of period T is characterized 
by a value of y in the range y to y + dy. This is given by 

P y (T,y)dy = \ T _71/2 y" 3/2 L(ky) dy . (3-13) 

We denote by F(T)dT the "population" of pulsars (sufficiently 
close and sufficiently luminous to be observed) with periods in the 
range T to T + dT. If f(y)dy is the number of pulsars which are 
observed to be characterized by values of y in the range y to y + dy, 
then 
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f(y) = f dT P y (T,y) F(T) . 


(3.14) 


On using equation (3*13) an d noting that N q , the total number of pulsars 
observed, is given by 


n o = f f(y)dy , 

O 


(3.15) 


we find that 


f (y) = N \ 1 k 1/2 y 3/2 L(ky) , 


(3.16) 


where 


X = J* dz z ^ 2 L(z) 


(3.17) 


A more detailed, generalized, derivation of equation (3-13) is made in 
Appendix B. In Appendix C we discuss the calculation of the function 
L(x) and its integrals. 


IV . Comparison with Observational Data 

Two sets of histograms have been constructed from observational data 

for comparison with the theoretical distribution given by equation ( 3 . 16 ), 

We adopt = I 3.5 km, which is close to the values computed by Cohen 

and Cameron ( 1971 ) for neutron stars in the mass range 0.3M„ ^ M„ < 2M . 

The theoretical histogram is taken from Appendix C; its midpoint is found 

from the relation x = 1 , i.e., log y = - log k, with k given 

pG 8 K peak 

by equation ( 3 .IO). Using the values of y computed from the data for 

each 7 ), the observed distribution is constructed for the same bins 

around y , that were used for the theoretical distribution, 
peak 
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To make Sample A as large as possible, we have included all 49 
pulsars for which pulse widths had been published at the time Paper I 
was written (Davies and Large 1970; Maran and Modali 1970) . To make 
this sample as uniform as possible, we have adopted for the width the 
value W quoted by Maran and Modali which, in some cases, is the "full 
width" and in other cases is the "equivalent width" . Sample B consists 
of the 13 pulsars for which Gunn and Ostriker (1970) list more careful 
measurements, and 5 additional pulsars from Maran and Modali (1970). 

These are listed, with our judgment of the best measure of W, in Table 

3 . 1 . 

Histograms formed from Sample A are compared with histograms generated 
from the theoretical distribution in Figure 3-l> for a range of values of 
T|. The best fit is found in the range 0.6 <, Tj ^ 0.7* The fact that the 
widths of the curves do not match too well is attributed to the intrinsic 
scatter in the measured pulse widths. Similar results are obtained from 
samples composed of those 25 pulsars with T > 0 .7 sec and of those 24 
pulsars with T < 0.7 sec. 

Histograms formed from Sample B are compared with the theoretical 
histogram in Figure 3-2, for a range of values of Tj. For this sample, 
the best fit seems to occur in the range 0.5 ^ 1| s 0,6. Although Sample 
B is smaller than Sample A, it gives a better fit to the theoretical 
curve. One may attribute the change in the value of Tj which gives the 
best fit, in going from Sample A to Sample B, to the fact that the total 
width is greater than the equivalent width. 

The above considerations suggest that T| = 0.6 ± 0.1. This uncer- 
tainty may reflect errors in the determination of Tj, or it may reflect 
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a real variation of Jj among pulsars and/or in time. We see from 
equation (3.4) that this leads to a braking index n = 2.2 ± 0.2. This 
quantity has been measured only for the Crab pulsar, for which various 
values have been proposed. Among the more recent publications, Richards 
et al . (197O) give n = 2.l6 ± 0.10; Duthie and Murdin (1971) give n = 
2.229 ± 0.003; and Roberts and Richards (1971) quote a range of values 
(for different time intervals) varying from n = 1.80 to n = 2.88. It 
appears that the value of n (and hence of T|) for the Crab pulsar may 
well vary, and that these changes in magnetospheric structure are related 
to irregular changes in period, but that the mean value of n for the 
Crab pulsar is not very different from that inferred for Samples A and B. 

V. Interpulse Phenomenon 

We may consider the number of pulsars which one might expect, on 
the basis of the above model, to show an interpulse. Geometrical consi- 
derations, discussed in detail in Appendix D, show that the probability 
P^(t) that an observer will see at least one pulse from a pulsar of 
period T is given by 

P X (T) = 3(R/R y ) 1/2 f sin 2 ® G( ©) d® = I.863 (R*/Ry) 1/2 , (3-l8) 

*0 

and that the probability P 2 (t) that the observer will receive radiation 
from both magnetic poles is given by 

P 2 (T) = | [l|l b (*/2)] 2 = 0.433 (VV • (3.19) 

Hence we find that the probability that an observed pulsar will show an 
interpulse is given by 
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p i “ P 2 /P l = °- 23 (VV 


1/2 


( 3 . 20 ) 


One may now estimate the expected number of pulsars, of a given 

observed sample, which should show an interpulse by evaluating P for 

each pulsar and summing these values. When this is done for Sample A, 

considering 7 ] in the range 0.70 2 71 2 0 . 50 , we find that the number 

of pulsars expected to show interpulses is in the range 0.8 to 1.7. 

This is to be compared with the fact that four pulsars from Sample A 

( PSR 0532 , 0904 , O95O, and I929) show interpulses. It may also be noted 

that three of these have short periods, consistent with the fact that 

P. «T - 1 ' /2 . 

1 


VI . Discussion 

We consider briefly some of the implications of the above results. 

6 1 

If, for definiteness, we adopt 7] = 0.6 and = 10°' cm, we find 
that, for T = 1 sec., R^ = 10^ cm whereas R^ = 10^'^ cm. Hence 
R^ is very much smaller than R^ . These considerations therefore indi- 
cate a pulsar magnetospheric structure very different from that envisaged 

by Gold ( 1969 ) ) Goldreich and Julian ( 1969 )? an d Gunn and Ostriker ( 1969 ). 

-is 7 Ip 

For the Crab pulsar, T = 10 sec., so R^. = 10 cm and R^ = 

8 2 -1 
10 * cm. On noting that the power extracted from the star is 2jtT N, 

and on equating this to the observationally indicated value of 10' 5 erg 


sec 1 (Shklovsky 1968), we find from equation ( 3 - 3 ) that = 10 J " L "^gauss, 

12 12 7 

which is considerably less than the values (10 - 10 ‘ gauss) which 

have been previously estimated. 

Scargle and Pacini (1971) have suggested that the "glitch" pheno- 
menon may be caused by events occurring in pulsar magnetospheres. This 


11.2 
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indicates a more complicated model of pulsar magnetospheres than has 
hitherto been contemplated. The results of the present study point in 
the same direction. In the next chapter we discuss a model magnetosphere 
which explains both the period-pulse-width distribution and braking 
index, and also offers a more detailed explanation of the glitch pheno- 
menon in terms of magnetospheric behavior. 
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TABLE 3.1 


The 1 8 pulsars making up Sample B, with their periods (t) and our judge- 
ments of the total pulse width (W) . Those marked in the last column 
(G and 0) make up the sample of 13 well determined pulse widths used by 
Gunn and Ostriker (1970). 


SAMPLE B 


Pulsar 

T( sec) 

W( sec) 

G & 

PSR 0329 

0.7145 

.050 

X 

NP 0527 

3.7455 

.090 

X 

NP 0532 

0.0331 

.006 

X 

PSR 0628 

1 .2444 

.140 


MP 0736 

0.3749 

.034 


CP 0808 

1.2922 

.070 

X 

PSR 0833 

0.0892 

.004 

X 

cp 0831 + 

1.2738 

.050 

X 

cp 0950 

0.2531 

.030 

X 

psr 1133 

1.1879 

.050 

X 

AP 1237 

1 . 3824 

.060 


PSR 1508 

0.7397 

.030 

X 

PSR 1749 

0.5626 

.020 

X 

PSR I 9 I 9 

1 . 3373 

0 

0 

X 

PSR I 929 

0.2265 

.016 


PSR 2016 

0.5580 

CO 

CM 

O 

X 

PSR 2045 

1.9616 

.100 

X 

AP 2303 

1.5759 

.030 
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FIGURE CAPTIONS 


Figure 3-1 • Distribution of Sample A of pulsars in the variable y, 

2— T) -2 

y = T 1 W , compared with histograms generated from 
the theoretical curve (broken line), for various values 
of the parameter 7 ]. The best fit occurs for 7 | in 
the range 0.60 ^ 7] ^ 0.70- The curves clearly do not 
fit for the value 7 ] = 1 . 

Figure 3*2. Distribution of Sample B of pulsars in the variable y, 

2-Ti -2 

y = T 1 W , compared with histograms generated from 
the theoretical curve (broken line), for various values 
of the parameter T|. The best fit occurs for 7 ] in 
the range 0.50 s Tj ^ 0 . 60 . 
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Figure 3.1 
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Chapter 4 


THE EFFECTS OF NON-RELATIVISTIC GAS IN THE MAGNETOSPHERE 

I . Introduction 

In most models of pulsar magnetospheres (Gold 1969; Goldreich and 
Julian 1969; Michel 1969 5 Ostriker and Gunn 1969; Gunn and Ostriker I97O) , 
it is assumed implicitly or explicitly that magnetic field lines con- 
tained within the "light cylinder" are closed whereas those which pene- 
trate the light cylinder are open. If the radiation from pulsars ori- 
ginates at the magnetic polar caps (those areas at the neutron-star sur- 
face penetrated by open field lines) and is beamed along the magnetic- 
field lines at these locations (Radhakrishnan and Cooke 1969; Komasaroff 
1970 ; Sturrock 1970 , 1971 a) there is a discrepancy between the inferred 
period-pulse-width distribution and the distribution derived from obser- 
vation (Gunn and Ostriker I97O; see Chapter 3) • This is illustrated, 
for Sample B of Chapter 3 > in Figure 4 . 1 ; the expected distribution for 
R Y= is shown as curve 1 . 

There is also a discrepancy concerning the braking index n (n = 

. o 

oxn/tt) ) : if the braking is due primarily to magnetic stresses, these 
models give n = 3; if the braking is due to a relativistic stellar 
wind, n = 1 (Michel 1969) . There is a great deal of scatter in the 
value of n derived from timing measurements of the Crab pulsar, but 
it is generally agreed that n > 1 always and n < 3 almost always; 
Boynton et al . (I972) have found that n = 2.4 ± 0 . 2 , and Rankin, Campbell, 
and Counselman (1971) have found that n^ 2.3 during quiet periods. 

II. Non-Relativistic Gas in Pulsar Magnetospheres 

The possibility that pulsar magnetospheres may contain substantial 

24 

amounts of non-relativistic gas (~ 10 g) was proposed by Scargle and 


33 



Pacini (1971). They introduced this hypothesis as an explanation of 
"glitches” (sudden changes in the periods of the Crab and the Vela pulsars) 
which, unlike the explanation in terms of "starquakes" (Ruderman 1969* 1972 ) 
would enable one to understand changes in the dispersion of the Crab 
radio pulses (Rankin and Roberts 1971 ) and in the optical appearance of 
the wisps in the Crab nebula (Scargle and Harlan 1970) which appear to be 
related to the September 1969 glitch in the Crab pulsar. The presence 
of a large amount of gas in pulsar magnetospheres, and fluctuations in 
this configuration, would offer an explanation of the prevalent "noise" 
in the residuals of timing data for the Crab pulsar (Boynton et al . 1972 ). 

The presence of such a large amount of cool gas in a pulsar magneto- 
sphere implies that gas is evaporated from the neutron star surface, pre- 
sumably by heating due to bombardment of the surface by high-energy 
particles. However, such a large mass of gas cannot be supported against 
gravity or contained against centrifugal force by the magnetic field. 

Hence appreciable mass accumulation can occur only in the region where 

these two forces are in approximate equilibrium, that is, at the "force- 

-2 

Rp 0 . Gas cannot collect at a 
radius appreciably less than R^. If if collects at a radius appreciably 
greater than R^ , centrifugal force will eventually fling out the 
collected gas, thereby opening the magnetic field lines to form two open 
flux tubes separated by a current sheet. This leads to the conclusion 
that Ry » Rpg. 

Ill . Braking Index 

-3 

In this model B « r for R^ ^ r ^ Rpg, whereas approximately 
B « r“ 2 for R^ B 

light cylinder is given approximately by 


* r < \> 


so that the magnetic field strength at the 


balance radius 


^B 


where 00 


^B - 


= GM„ 
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(4.1) 


b l - B A 3 hi 1 > 
where B„ is the magnetic field strength at the surface of the star. In 

*7v 

previous calculations of the magnetic torque exerted on a pulsar, it has 
been assumed that the toroidal and radial components of the magnetic 
field are comparable at the light cylinder. It seems reasonable to 
retain this assumption in the revised model, since if relativistic gas 
is flowing substantially radially outward at the light cylinder, and is 
also moving along magnetic field lines, the spiral angle must be approxi- 
mately U 5 0 at the light cylinder. (Similar arguments concerning the solar 
wind lead to the accepted formulas for the magnetic torque on the sun, 
except that is replaced by the radius at which the solar wind speed 

is equal to the Alfven speed (Dicke 1964; Weber and Davis 1967 ) . Hence, 
using equation (2.6) of Sturrock ( 1971 a), 

» - 5 V s < ( 4 - 2 > 


one obtains the following expression for the magnetic torque: 




(4. 3) 


7/3 

Since Nad) , we see that n = 7 / 3 > in good agreement with the value 
n « 2.4 ± 0.2 derived from analysis of timing data of the Crab pulsar 
(Boynton et al. 1972). 


IV. Period-Pulse-Width Distribution 

We may compute the period-pulse-width relation to be expected in 

the present model. The values should cluster about the relation W/T = 

3 1/2 

2 l|l^>/ 2 jt, where ^=2 ’ This leads to the relation 
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W m 10 1 - 15 M^" 1/6 R* 1/2 T 273 




This relation is shown as curve 2 in Figure Ij-.l, adopting NL = M and 

O 

6 l 

R v = 10° ' cm. The overall agreement is good. 

Following the lines of Chapter 3 and the generalization of Appendix 

2 1/3 

B [with T| = 2/3 and R^ = (( M^/kn ) , we may estimate the scatter 

to be expected in the period-pulse-width distribution due to random 

orientations of the line of sight and of the magnetic and rotation axes. 

The resulting theoretical distribution of pulsars in the variable 
2/3 2 

y = (T ~yw) is compared with the observationally determined distribu- 
tion, for a range of neutron-star models (Cohen and Cameron 197l)> in 
Figure 4.2. The agreement is good for neutron-star masses in the range 
0-5 M q <; ^ M q . 

V. Glitches 

2 

During active periods, the moment of inertia I (g cm ) of the 

star plus magnetosphere will be varying. At all times, in addition to 

the magnetic torque, there may be a loss of angular momentum at a rate 
2 ”2 

L (g cm sec ) due to the escape of particles from the system. These 
effects lead to the following expression for n: 


, a) I cd 2 I 1 , . dL N 

n-n =(n - 2) — — iiL-co — ), 

o v o ' • I *2 1 ,vr ° dco ' ’ 

CD CD CDi 


(M) 


where n = t; 4— = ■?. A glitch is now interpreted as an instability of 

o N dcD 3 

the gas-magnetosphere configuration which leads to the ejection of some 
of the gas from the magnetosphere and/or the dumping of some of the gas 
to the surface of the star. There is therefore a change 51 in the 
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moment of inertia and a change 6H in the angular momentum of the star- 


magnetosphere configuration. This leads to a change in the angular velo- 
city given by 

JS03 = _ 61 + 6H 
cn I cnl ’ 

For the process envisaged, both 61 and 6H are negative, which allows 
6co to be either positive or negative. 

After a glitch, the moment of inertia will return slowly to its 
original value. If the rate of evaporation from the neutron-star sur- 
face is approximately constant, and if the rate of mass loss from the 
collected gas is proportional to the mass of gas, the moment of inertia 
of the magnetosphere will relax according to the law 

Al(t) = AI[1 ~ e _t/TG ] . (4.7) 

Noting that immediately prior to a glitch the moment of inertia of the 
magnetosphere is Al> we find that (despite the fact that o>T G » 6co) 
the post-glitch angular velocity m^ G may be related to the extrapolated 
"no-glitch" angular velocity <n^ G by 

= 1 + M + M e‘ t/T G . (4.8) 

col I 



The post-glitch braking index np G is related to the steady braking 
index n by 


n 


PG 


n 


o£ AI J_ -t/T G 

• 2 1 2 
0) r G 


(M) 
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Equations (4.8) and (4.9) are in fair agreement with the post-glitch be- 
havior of the Crab and Vela pulsars. 

The September 1969 glitch of the Crab pulsar comprised a sudden 
increase in rotation rate by three parts in 10^. The subsequent relaxa- 
tion of the angular frequency agrees approximately with equation (3.8) 
with 6H = 0 and t g approximately two months (Roberts and Richards 
1971 ; see especially Figure 2). The observed value of 6aV’co 
(10 ; Papaliolios, Carleton and Horowitz 1970) leads to a 

similar estimate for t^,. Equation (4.9) gives an approximate represen- 
tation of the behavior of the braking index which was close to 2 prior 
to the glitch, jumped to nearly 3 after the glitch, and subsequently 
decayed to ~ 2.5 in about six months (Roberts and Richards I 97 I ; Table l) . 

The change fica/cn » 10 requires the dumping of 10^ * g of gas 

from the force-balance region to the neutron-star surface. The potential 

or kinetic energy of such a mass in orbit at the force-balance radius is 

4o 3 

of order 10 erg, comparable with the energy required to produce the 

observed wisp disturbances (Scargle and Harlan I970) . The inferred 

post-glitch time scale of two months implies an evaporation rate of 

1<5 r -1 35 4 -1 > 

10 g sec which requires a power of 10 erg sec , about 0.270 

of the power developed in the Crab pulsar. 

Preliminary calculations indicate that this model of a glitch is 

consistent with published data (Radhakrishnan and Manchester 1969 ; 

Reichley and Downs 1969 * 1971) concerning glitches in the Vela pulsar, 

but more definite conclusions must await the publication of more detailed 

observational data. 
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VI . Period-Age Distribution and Magnetic Fields 

The period-age distribution of pulsars (Sturrock 1971b) may be 
reconsidered in terms of the present model. The results are shown in 
Figure 4 . 3 . It is seen that the "turn-off" condition associated with 
the cessation of pair-production in the electron polar zone is unchanged 
and in good agreement with observations. We still find that the Crab 
pulsar is the only one spinning rapidly enough for pair-production to 
occur at the ion polar zone, which offers an explanation of the fact that 
this is the only pulsar exhibiting a precursor. The most important 
Ghange is that the present model leads to estimates of smaller by 

a factor of 10 than those based on earlier models. This change goes 
some way towards resolving the discrepancy noted by Gunn (1971) con- 
cerning the relative magnetic field strengths of neutron stars and 

11 3 

white dwarfs. For the Crab pulsar, we find 10 gauss, which is 

12 12 6 

to be compared with earlier estimates of 10 to 10 * gauss. 


39 



FIGURE CAPTIONS 


Figure 4.1. Period-pulse-width distribution for Sample B of Chapter 3* 
Curves 1 and 2 are theoretical curves expected on the basis 
of the polar-cap model for Ry = R^ and Ry = , res- 

pectively. 


Figure 4.2. Distribution of Sample B (Chapter 3) of pulsars in the variable 

2/3 2 

y, y = [T /W] , compared with histograms generated from 
theory (broken lines) for Ry = and various values of the 

neutron star mass (Cohen and Cameron 1971) • 


Figure 4.3. Distribution of pulsars in period and "age" [ t=T/( dT/dt) ] 
compared with theory. Q, S-type pulse structure; ZS. , 
C-type structure; □ , D -type structure, in the notation 
of Taylor and Huguenin (I971). 
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Figure 4.3 






Chapter 5 


DETAILED MODEL OF THE COLLECTION OF GAS 

I . Introduction 

A model of pulsars in which the inner magnet ospheric structure is 
as discussed in Chapter 4 seems to resolve two major discrepancies between 

for which 

data. It is therefore desirable to analyze this model in more detail. 

More complete knowledge of the magnetospheric structure would make it 
possible to improve the calculation of the theoretical period-pulse-width 
distribution and to estimate the current distribution in the polar caps 
which would then provide a basis for calculation of pulse shapes. Calcu- 
lation of the dynamic properties of the magnetosphere would provide infor- 
mation which could be compared with irregularities in timing data such as 
the persistent "noise” and "glitches." We shall find that this model also 
suggests an explanation of the "drifting sub-pulse" phenomenon (Hewish 
1970) which is a puzzling feature of pulsars. 

If a pulsar magnetosphere is initially free of gas, so that is 

approximately and gas evaporates from the surface of the star due 

to the bombardment of the surface of the star by relativistic particles, 
we can determine the locations at which gas may initially collect. This 
is a question of microscopic stability, in which we may for simplicity 
consider a single particle constrained to move on a corotating magnetic 
field line under the influence of gravitational and centrifugal forces. 
This problem is set up in Section II. If there is a progressive accumu- 
lation of gas, there will at some stage be sufficient gas in the trapping 
zones that the configuration will become unstable. This question of 


the properties of earlier models ( 


V 


and observational 
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macroscopic stability is much more complex and will not be discussed in 


detail in this thesis. Instead, we discuss in simple physical terms 
the motion of collected gas in the limit that the mass becomes very large. 
These physical arguments show that matter collected substantially beyond 
the force-balance radius will be flung outward, opening the magnetic field 
lines to which it is tied. 

In Sections III and IV, we discuss the aligned and orthogonal cases 
in detail. In each case, there are small collection zones which are 
both microscopically and macroscopically stable. By determining where 
macroscopic stability changes to instability, we determine the separa- 
trix between closed and open field lines, and hence determine the sizes and 
shapes of the polar caps. With this information, it is possible to re- 
examine the period-pulse-width distribution and to compare it with observa- 
tional data; this is done in Section V. The interpulse phenomenon is exa- 
mined in Section VI, and further topics are discussed briefly in Section VII 

II . The Model 

To investigate the trapping of gas and its effect on the structure 
of the magnetic field, we find the locus of stable points for single 
particles constrained to move along field lines. The magnetic field 
(surface strength B # ) is assumed to be dipolar from the neutron star 
to the light cylinder. Thus we have the problem of "a bead on a wire”; 
the bead is subject to gravity, and the wire rotates with a constant 
angular velocity. The stable points of the motion will tell us where 
cool gas will begin to collect in the magnetosphere. 

In a uniformly rotating reference frame, the Lagrangian for a 
particle of mass m is given (Landau and Lifshitz i960) by 
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(5.2.1) 


_ 1 2 -4 -* "♦ 1 2 „ /-4\ 

L = — mv + mv • a) x r + - m(<x> x r) - U(r) , 

where ufr) is the potential energy of the particle in an inertial frame. 
The motion is one-dimensional so the term mv • a) x r makes no contri- 
bution to the Euler-Lagrange equations and the particle moves with effec- 
tive potential energy 


V = - | (£D x r) 2 - CM* m/r . (5-2.2) 

The potential energy U(r) is the gravitational potential energy due 
to a spherically symmetric central body of mass M*. To evaluate V, 
we define coordinates in the rotating frame such that the magnetic dipole 
axis is in the z-direction and the angular velocity vector is in the 
x-z plane with polar angle ®. If ( r, 0,<j>) are the usual spherical 
polar coordinates, 


v = -t 


2 2 
<cjo r 


[si 


■ 2 . 
sm 0 


sin 2 <|>+ |sin0 cos<j) cos® 


COS0 



caim 

Tv 


(5-2.3) 


In dimensionless form, all lengths are scaled to the force-balance 
radius R^, all energies are given in terms of a typical energy E q , 
and all forces in terms of a typical force F q where 

Rp B = (cm*cd' 2 ) 1/3 , (5.2.4) 


E 

o 


2 t>2 

m 0) R^b 


y 


(5-2.5) 


F 

o 


m cn 2 RpB . 


( 5 . 2 . 6 ) 
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The dimensionless quantities are defined by 


- °*TB 


, r = BR,^ ; 
’ o F FB 


(5-2.7) 


v = ye . E = eE ; 
o 3 o 


( 5 - 2 . 8 ) 


F = f F . 
o 


(5-2.9) 


The resulting expression for the effective potential is 


Y = - ~~ jsin 2 0 sin' 


in 2 (j) + |sin 0 cos(j) cos® - cos 0 sin®j 2 + 2 OL .( 5 . 2 . 10 ) 


Furthermore, we may label a particle by the magnetic field line to 
which it is tied. A dipole field line is given by 


<j) = constant , r 1 sin 2 0 = r ‘ L , 


-1 


(5-2.11) 


where r Q is the greatest extent of the field line. In dimensionless 
form, 


2 

(|) = constant , p sin 0 = a , 


( 5 . 2 . 12 ) 


so that (f 3 ,<|)) identifies a field line and 0 locates a particle along 
it. For a particle so labeled, Eq . ( 5 . 2.10) becomes 

12 4 1 2 p / \p _o — ^ / 

Y = - ^ P sin Bjsin^© sin (j)+ Isin 0 cos(j) cos©-cos 0 sin®n+ 2 f 3 _ ' 5 sin 0/ . 

(5.2.13) 

The rate of change of Y along a field line is given by 


dY/ds = n • VY , 


(5-2.14) 


where n is the unit vector parallel to the magnetic field; 


? — I/O 

n = (1 + 3 cos 0) (2 cos0, sinQ, 0) 


(5.2.15) 


From Eq. (5.2.10 ) , 


R pB (dY/ds) = a(l+3 co s 2 e) 1//2 |-2cos0^sin 2 0sin 2 (|) + |sin0cos(j)cos®- cos0sin®j 2 j 

- sin0 ^sinOcos 0sin 2 <j> + ^sin0cos<j)cos® - cos0sin®j 

X |cos0cos(j>cos® + sin0sin®jj + 2cos0 a . (5.2.16) 

Along a field line labeled by (|3,(|)), this is 

R^dVds) = psin 2 0^1+3cos 2 ej 1/,2 |-2cos0^sin 2 0sin 2 (jH- |sin0cos(|>cos®- cos0sin®j 2 J 

- sin0^sin0cos0sin 2 (j) + |sin0cos<j)cos® - cos0sin®| 

X |cos0cos<j>cos® + sin0sin®jj + 2(3 ^cosOsin ^q| , (5.2.17) 

Another quantity of interest is the energy of a particle tied to a 
field line, as measured in an inertial frame (at rest with respect to the 
origin of the rotating frame). Neglecting the kinetic energy due to small 
velocities in the rotating frame, 


E = | (u> x r) 2 - GM # m/r . 


( 5 - 2 . 18 ) 


In dimensionless form, this is 


e = ^ o! 2 ^sin 2 0sin(J) + | sin0cos({>cos® - cos0sin®j ? - 2C£ . ( 5 . 2.19) 

Along a field line labeled by (p,<b). 


g = i p 2 sin^0^sin 2 0sin 2 (j) + |sin0cos(j)cos® - cos0sin®j 2 - 2(3 ^sin 


( 5 - 2 . 20 ) 
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Finally, we need the force on a particle in the rotating frame. 


From Eq. (5.2.2), 


? = -nS x (ffl X r) - Cai„mr/r^ 


(5-2.21) 


In dimensionless form, this has components 


f^ = oi|sin^8sin^(|> + |sin0cos((>cos® - cos0sin®j^ - a , (5.2.22a) 

= o:|sin0cos0sin^(j) + |sin0cos<|>cos® - cos0sin® j|cos0cos(|)cos®+ sin0sin®jj , 

(5.2.22b) 


f^ = a|sin0sin<j>cos<|) - sin(j)cos® |sin0cos<j)cos® - cos0sin®jj . (5.2.22c) 


For a given field 
These results may 
Considerable 
formations 


line, we may simply substitute according to Eq. (5.2.12). 
be confirmed by differentiation of Eq. (5-2.10). 
simplification of our work occurs because the two trans- 


<|> — ► 2jt — 

and 

<j) -* Jt - 
0 -* ir - 

leave f unchanged (see Eq. [5.2.13])* 
only the intervals 0 ^ 0 ^ jt/2 and 0 
will be used throughout: 


♦ (5.2.23) 

(5.2.24) 

Thus we need explicitly consider 
£ £ it/2. The following notation 



X = 

2 

sin 0 , 

(5.2.25) 

n = 

2a 

cos <]) , 

( 5 . 2 . 26 ) 

X = 

2 

cos ® . 

(5.2.27) 
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Note that a single value of the coordinate corresponds in general to 
two values of the coordinate 0 in the complete range 0 ^ 0 ^ it. 

The requirement for microscopic equilibrium at a given point is 

dy/ds = 0 . ( 5 . 2 . 28 ) 

The condition for microscopic stability at that point is that, in 
addition, 

d 2 Y/ds 2 > 0 . ( 5 . 2 . 29 ) 

However, it is not necessary to use condition ( 5 . 2 . 29 ) to decide which 
equilibrium points are stable. Since ¥-*-«> as r -* 0, there must 
be more than one solution of Eq. (5.2.28) if Eq. (5.2.29) is to hold. 
Specifically, for a given (p,(f»), a stable point 0 s (f3,<j>) is a solution 
of Eq. (5.2.28) which lies between two other solutions of this equation 
which correspond to values of y greater than f( 0 g ) . 

This rule has a simple physical interpretation. On a field line 
which originates at low magnetic latitude, there is only one equilibrium 
point and this is an unstable point because the inward gravitational force 
dominates the centrifugal force. On magnetic field lines which originate 
at large magnetic latitudes, there are three equilibrium points: the 

outermost point is in a region where the outward centrifugal force domi- 
nates, so it is stable; the other two points are where the components of 
gravitational and centrifugal forces along the magnetic field cancel out, 
and each of these is an unstable point . Hence if one knows the number of 
equilibrium points on a given field line, one can immediately determine 
whether each point is stable or unstable. Application of these rules will 
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show that, for any alignment between the magnetic and rotation axes, 
stable collection points exist on field lines which extend out approxi- 
mately to or beyond. 

We now consider briefly the question of macroscopic stability, i.e., 
the stability of a mass sufficiently large that the concomitant dynamical 
and gravitational stresses are comparable with the magnetic stresses. A 
detailed treatment of this topic represents a complex problem of magneto- 
hydrodynamics which must be deferred for further study. In discussing an 
onset of instability, one of the key questions is to determine whether the 
onset is explosive or non-explosive (Sturrock 19 66) . This involves know- 
ledge of the non-linear processes involved, but the question may often 
be decided without making a complete non-linear dynamical analysis. 

What happens as mass begins to accumulate at microscopically stable 
points? There are three regions we must consider: 

(l) Material collects where the radial force is inward (and the 
energy E is negative ) . This can only occur inside and close to the 
force-balance radius. Were mass in these places decoupled from the magne- 
tic field, it would orbit with angular velocity somewhat higher than that 
of corotation. This implies that, as mass builds up, it will stretch 
the magnetic field lines. The magnetic stress will enforce corotation 
of masses smaller than some critical value, but larger masses will cause 
the magnetic field to wind up progressively. Since this would lead to a 
progressive increase in the transverse pressure, the configuration is 
likely to become MHD unstable to an instability of the Rayleigh-Taylor 
type (Schmidt 1966) . A competing instability, in this situation, would 
be the tearing-mode instability (Furth, Killeen and Rosenbluth 1963) 
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which would lead to field-line reconnection. The MHD instability is 
likely to be explosive and may cause some material to be ejected and 
some to be dumped to the surface of the star, but the resistive insta- 
bility is probably non-explosive. 

(2) Material collects where the radial force is outward but the 
energy is negative. This can only occur outside and close to the force- 
balance radius. This case is similar to Case (l) except that the magne- 
tic field will tend to wind up in the opposite sense since material 
would tend to orbit with angular velocities lower than that of coro- 
tation . 

(3) Material collects where the radial force is outward and the 
energy is positive. This is the case for all collection points beyond 
the force-balance radius, except for the small region of Case 2. As 
the mass builds up, it stretches the field lines outward as it attempts 
to enter an open orbit. The magnetic stress will enforce corotation 

of masses smaller than a critical value, but larger masses will enter 
open orbits, carrying the attached magnetic field lines with them. This 
leads to an open magnetic field configuration. The transition appears 
to have the character of an explosive instability. 

The implication of the above considerations is as follows. If 
the magnetosphere is initially free of non-relativistic gas and evapora- 
tion occurs, gas will collect on field lines which extend to approxi- 
mately Rpg or beyond. When a sufficiently large mass has accumulated, 
the magnetic field will be distorted. On field lines denoted by Case ( 3 ) 
there will be a change of configuration (possibly explosive), leading to 
the ejection of gas and the conversion of closed field lines into open 
field lines. On field lines denoted by Case (2) the change in configura- 
tion would not initially be a drastic one, although the accumulation of a 
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sufficiently large mass of gas might eventually lead to a disruption 
of the magnetic field, possibly converting closed field lines to open 
field lines. On the field lines denoted by Case (l) the initial change 
in configuration will not be drastic, although a sufficiently large mass 
may lead to a large distortion or disruption of the magnetic field con- 
figuration. It seems quite possible that both the glitches and the 
"noise" observed in the Crab pulsar (Boynton et al . 1972 ) are the result 
of processes such as the above occurring in the pulsar magnetosphere. 

We also note the possibility that trapped material in the force- 
balance region may be coupled to magnetic-field lines which have become 
disconnected from the main magnetic-field pattern of the neutron star, 
due to differential rotation leading to field-line reconnection as men- 
tioned above. If this occurs, and if the trapped material has an azi- 
muthally nonuniform distribution, the magnetosphere will be subject to 
a perturbation drifting azimuthally with respect to the neutron star. 
Such a perturbation would produce an azimuthally drifting perturbation 
of the current distribution at the polar caps. Since the radio emission 
is due to these currents, the radiation would display a periodic modu- 
lation similar to the observed "drifting subpulse" phenomenon. 

In what follows, we are concerned primarily with the separatrix 
between closed and open field lines. To this end, it is convenient to 
introduce the following criteria, corresponding to the inner and outer 
boundaries of Case 3: 

Criterion 1. Radial force zero or outward. 

Criterion 2. Energy zero or positive. 

The preceding considerations indicate that Criterion 1 gives an inner 
limit to the separatrix, and Criterion 2 gives an outer limit. In 
analyzing the polar-cap configuration and the period-pulse-width distri- 
bution, we shall give results for both criteria, but it will be seen 
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that it makes little difference which one is adopted. 

As a simple application of the preceding rules, we discuss briefly 
the useful case of a field line at <j) = ft/ 2 or <|> = 3it/2 ( T} = 0), and 
show that there exist stable points for all ©. In this case Eqs . 
( 5 . 2 . 17 ), (5-2.20), and (5-2. 22a) reduce to 

R FB ( d ^ /ds ) = PX(^-3X)“ 1 / 2 (1-X) 1 / 2 G(X) , G( X ) = -3XX " 2(l-X) + 2p"V 3 

(5-2.30) 

e = \ P 2 X 2 H(x) , H( x ) = XX + (1-X) " 2p“V 3 (5-2.31) 

f r = PXK(X) , K(x) = XX + (1-X) " P"V 3 • (5.2.32) 

From the first equation it follows that \ = 1 is an equilibrium point 

for all p and It is a stable point if there is a single root of 

G(x) =0 in the interval 0 < X < !• As 

G'(x) = -3X - 6 p"V 4 , (5.2.33) 

there are no extrema of G in the interval. Thus there is only one 
root of G(x) = 0; this appears, initially at x = 1, when G(l) ^ 0. 

Therefore the point (x = 1, p, T) = 0) is a stable point if 

(3 > P S (H = 0,X) = [2/(2fX)] 1/3 . (5.2.34) 

Next we consider the energy and radial force at these collection 
points. It is found that e ^ 0 if 

p £ p e (n = 0,x) = 2 1/3 « 1.260 (5.2.35) 

and that f S 0 if 
r 
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(5-2.36) 


p * p f (n =. o, x) = i . 

Thus stable points occur at x = 1 whenever (5.2.34) is satisfied, 
the energy is zero or positive at these points when ( 5 - 2 . 35 ) holds, and 
the radial force is zero or outward if ( 5 . 2 . 36 ) is met. Hence we have 
demonstrated that stable collection is always possible for those field 
lines at <j) = jt /2 and <j) = 3 ^/ 2 , and that whichever criterion is 
adopted for the opening of field lines, it occurs at a distance which 
is independent of 0. If criterion (l) is accepted, then, apart from 
the case © = tc/2, there is a region in which particles may collect and 
remain without opening the magnetic field pattern. If criterion (2) is 
adopted, such field lines exist for all ®. For either criterion there 
is a much larger region in which material can collect stably and will 
drag open the magnetic field. It will become clear in the following sec- 
tions that there is a finite range of <|) around <j) = n/2 and around 
(|) = 3jt/2 in which material can collect and be contained by the field. 

Ill . The Aligned Case 

For 0=0, Eqs. (5.2.17), (5.2.20), and (5.2.22a) reduce to 

Rpgdf/ds = Px (4-3x) _ 1 / 2 (l-x) 1 / 2 (-3x + 2ffV 3 ) (5.3.1) 

e = (1/2)P 2 X 2 (X - 2p"V 3 ) , (5-3.2) 

and 

f r = 6 x(X"P"V 3 ) , (5.3.3) 

respectively. From Eq. (5«3.l), X = 1 is an equilibrium point for all 
0 (all (f) are equivalent) . The other root of dy/ds = 0 on 0 < X < 1 
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is 


X = (2/3 P 3 ) lA (5.3.4) 

and since x < 1* this exists only if 

P > P s (® = 0) = (2/3) 1/3 « 0.874 . (5-3.5) 

From Eq. (5-3.2) the region of zero or positive energy for a given (3 is 
found to be 

X £ ( 2/ p 3 ) 1 A (5.3.6) 

and only those field lines which satisfy 

P * P e (® = 0) = 2 1/3 » 1.260 (5. 3- 7) 

have such a section. The radial force is zero or positive for a given 

P if 

x * P _3A • (5-3.8) 

Only field lines with 

P * P f (® = 0) = 1 (5.3.9) 

can satisfy Eq. ( 5 . 3 . 8 ). 

We have drawn this picture: All field lines outside of p = 0.874 

have microscopically stable points at x = 1 ( 9 = rt/2) . The radial force 
is outward for those stable particles outside of (3=1, and the energy 
is positive for those outside of p = 1.260. Thus for either criterion 
for the opening of field lines there is a region of macroscopic stability. 
In Figure 5 .1 we show the various regions of collection. Figure 5-2 
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is a comparison, approximately to scale, of the magnetospheric structure 
given by this model and that of Goldreich and Julian ( 1969 ) . 

It is useful to know the depth D(p)E Q of the effective potential 
wells around collection points. We find 


D(p) = (l/2){3 2 + ff 1 - (V3)(3/2) 1/ V lA ? P * (3 S (® = 0)] 


D(p) = 0 , 3 < p s (® = 0) . 


(5.3.10) 


D(1.26o) 0.195? so the potential well in the trapping region can be 

almost 20 percent of E q . Since the work W required to lift a particle 
from the surface to the collection regions is approximately W = GM^m/R^ 

= (Rp /R^)E , this depth is around 2 percent of W. However, one need 
only go as far out as p M 5 to reach D(p)E Q « W, a condition which 
makes collection of particles (and therefore opening of field lines) very 
likely. 


IV . The Orthogonal Case 

We now consider the case ® = jt/2. There is substantial evidence, 
including the swing of the polarization vector during a pulse of the Vela 
and other pulsars (Radhakrishnan et al . 1969 ; Manchester 1971b) and the 
existence of interpulses for six pulsars (PSRS 0531, O 923 ? 0904, 0950? 
1055, ’and 1929; Backer et al . 1973)? that the dipole axis is highly 
inclined to the rotation axis for these pulsars. When ©= jt/2, we 
have 

R FB dY/ds = Px(l"X) 1/2 (4-3 x)" 1/2 Q(x) ? (5- 4. la) 
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Q(x) = 3T1X - 2 + 2(3 V 3 ; 


(5.4.1b) 


e = | P 2 X 2 R(X) , 


(5.4.2a) 


R (x) = -T1X+ 1 - 2p“ 3 X _3 


(5.4.2b) 


f r = PXS(X) ; 


(5.4.3a) 


S(x) = -T 1 X + 1 - P“V 3 . (5.4.3b) 

Inspection of Eq. (5.4.1a) shows that X = 1 is a zero of dY/ds for 
all azimuthal angles (characterized by Tj) . Finding the other roots of 
dy/ds =0 is equivalent to solving 

Q( x ) = 0 , (5-4.4) 


which, in turn, is equivalent to solving 

x 3 Q(x) = 0 . ( 5 . 4 . 5 ) 

Although Eq. (5*4.5) is quartic in x and thus (in principle) soluble 
analytically, it is possible to obtain all the essential features of the 
solution without doing so. We only need decide which portions of (Tl,p)- 
space correspond to zero, one, and two roots of Eq . (5*4.4) in the range 
0 < \ < 1. Note that Q r (x) is zero in only one place, at 


and that q"(x) 


*o = 2 


1/4 -3/4 -1/4 

P T] 


(5-4.6) 


is always positive, so 


Q(Xq) 


is the absolute minimum 
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of Q(x) • If Q(Xq) I s negative there are two real roots of Eq. (5-4.4); 
we must discover how many, if any, of these lie in the desired interval. 
This may be done and the result is Figure 5 . 3 a, which shows the division 
of ( T], p) -space we sought. 

We have found the following: For any T| there is a p sufficiently 

large that there exists a microscopically stable point x„ at that 

azimuth. When J| < 1/2, the stable point is x„ = I» this requires a 
1 s 

l/o 

field line with p > [2/(2-3l))] • For Tj > 2/3; the stable point 

occurs at x <1; here we need p > If 1/2 < T| < 2/3, there 

is initially (p > 2r J T|) a stable point at some Xg < 1; as P in- 
creases the stable point moves toward x = 1; an d reaches x = 1 

s s 

I/O 

when p ^ [2/(2~3Tl)] • [When a stable point occurs at ^ < 1, there 

1/2 

are physically two stable points, one at 0 = arc sin(x ), and the 

s 

1/2 

other at 0 = it - arc sin(y g ).] 

Next we investigate the energy and radial force. From Eqs. (5.4. 2a) 
and (5.4.3a) it follows that the sign of g is determined by the sign 
of R(x); the sign of f by the sign of S(x) . The procedures out- 
lined above, applied to R and S, enable us to construct Figures 5-3 b 
and 5 . 3 c, which show the signs of these quantities in various regions of 
(TUP) -space. 

Comparison of Figure 5-3 a with Figures 5 .3b and 5-3 c is sufficient 

to determine the character of the stable points (y^ = l) for T] < 1 / 2 . 

Figure 5 . 3 c shows that f > 0 for all "f|. Figure 5-3 b shows that if 
1/3 

p < [2/(l-T])] J the energy is negative at the collection points. We 
also find that, for 1/2 < T] < 2 / 3 ; the stable points which occur at 
X g = 1 always have g > 0 and f > 0. 
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For Tl > 1/2 these results are not sufficient to determine the 
signs of e and f^ for those stable points located at x g < 1 . In 
order to do so we require more detailed knowledge of the roots of 0.(x) > 
R(x), and S(x) . The asymptotic form of the larger and smaller roots may 
be found by considering the limits p -• eo and X~* 0, respectively. 

In an obvious notation, 


xj - (2/3)(l/T) , 


(5.4.7a) 


Xq - (1/P) 


(5.4.7b) 


+ 



(1/Tl) , 


(5.4.8a) 


X^-2 1/3 (l/p) ; 


(5.4.8b) 


Xg - (Vn) t (5.4.9a) 

Xg - (i/p) • ( 5 . 4 . 9 b) 

The smallest possible real zero of Q(x) is subject to Q(x o ) = 0; 

this is 

Xq^ = 2 2/3 (l/p) = (1/2) (l/Tl) . (5-4.10) 

Similarly, the smallest real zeros of R and S are 

x ; ^ = 2 (i/p) = ( 3/4) ( i/*n) (5-4.li) 

and 

y 0-) = 2 2/3 (1/p) = ( 3/4 )(i/ti) . (5.4.12) 
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is a root of R(x) = 0, 


In addition, it may be verified that 



and that it lies on the branch of smaller roots. 


From these facts it follows that for all Tl. the region f s 0 

- r 

completely contains the region e S 0, and, furthermore, it contains 
the entire locus of stable points . In addition, only for T| < 1/2 do 
stable points occur where e < 0. These features are illustrated in 
Figure 5-4; numerical solutions for the zeros of Q, R, and S have 
been used to draw the curves . 


In the orthogonal case there are microscopically stable points at 
all azimuths. At none of these is the radial force inward. Stable 
particles (at 0 = it/ 2) in the intervals jt/4 ^ <j> ^ 3^/4 and 5^/4 <, 

2 Wo 

<fl ^ 7it/4 have negative energy if they are inside a. = (3 = [2/(2~3cos (J>)] , 

and will be macroscopically stable if Criterion 2 determines the opening 
of field lines. Collection of substantial amounts of material at all 


other microscopically stable points will lead to an open magnetic field 
pattern beyond those points. 


V. Polar Caps and Period-Pulse-Width Distribution 

The polar caps are those areas at the magnetic poles which are 
bounded by the intersections of the outermost closed field lines with 
the surface of the neutron star. The symmetry of 7 under the trans- 
formations of Eqs . (5.2.27) implies that the outermost closed field 
lines, P c (<}>)i obey 

P c (<j>) = P c (2it-<j>) (5-5-1) 

and 
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p c ( 4>) = e c (l H>) • 


( 5 . 5 . 2 ) 


This means that the shape of a cap is unchanged under rotations of l80° 
about the projections of the x- and y-axes onto the planes z = and 
z = -R*. 

If radiation is produced near the neutron star surface and is beamed 
along open field lines, the half-angle of the emission cone is approximately 

♦(♦) = ( 3 / 2 ) e r (4») (5.5.3) 

where 

e r («t>) « sine r ((}>) = (VR^) 172 f3 c ‘ 1/2 (<|>) . ( 5 . 5 .^) 

Since we will find that the polar caps and emission cones are roughly 
elliptical, we will speak of their semi-major (<j) = it/2, 3^/2) and semi- 
minor (<j) a 0, it) axes. 

In Chapter 4, it was assumed that the last closed field line at a 

given azimuth was the one which grazed a cylinder of radius R centered 

FB 

1/2 

on the rotation axis. This leads to semi-major axis ^ = ( 3/2) (R^/R^g) 
and semi-minor axis ^ = \j;^G(©). The definition and properties of the 
function G(®) are given in Appendix A. We introduce a reduced measure 
of the angular size of polar caps and emission cones, 

h(<}>) = P c _1/2 (<f>) y (5.5.5) 

so that the model of Chapter 4 has h(^/ 2) = 1 and k( 0) = G( ©) (see 
Appendix A) . 

Choosing alternatively Criterion (l) or (2) for the opening of field 
lines, we find from the example in Section II that the semi-major axis is 


62 



(*/ 2 ) = 1 


(5.5.6a) 


or 

k (2) (jt/2) = ( 1/2) 1/6 « 0.891 , (5.5.6b) 

respectively. We have no general expression for the semi-minor axis, so 
we discuss the two special cases. The aligned case involves a circular 
polar cap, with the size given by Eqs . ( 5 . 5 . 6 ). In the orthogonal case 
we determine the shape of the polar cap as follows: If T| > 1/2, collec- 

tion always occurs with positive e and f , so the first microscopically 
stable point determines the opening of field lines: 

P c (<j>) = 2 573 cos 2 <j> , <j> <; «A . (5.5-7) 

For 71 < 1/2 the collection occurs at x = 1 with initially negative 

s 


energy and positive radial force. Thus 

p ( X ) = [2/(2-3cos 2 <j >)] 1/3 , <t> ;> *A ; (5.5.8a) 

p( 2 ) = 2 l/3 sin _ 273 <J> , <j> ;» jt/4 . ( 5 . 5 . 8 b) 

Therefore, in the orthogonal case, the polar cap has the shape 

k(<| > ) = (l/ 2 ) 5 / 6 (l/cos<J>) , <{> <; */4 ; (5-5-9 a ) 

K (l) ((j>) = ( 1 / 2 ) 1 / 6 ( 2-3 cos 2 *) 176 , <j) 2 */4 ; ( 5 . 5 . 9 b) 

K^(<j>) = ( 1 / 2) 176 sin 1 / 3 (j> , <{> 2 rt/4 . (5-5-9 c ) 


We illustrate in Figure 5-5 the shapes' and relative sizes of the polar 
caps in the aligned and orthogonal cases. It is interesting to note 
that in the latter case the simplicity of the model has caused two sides 
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of the cap to be straight lines. Numerical examination of intermediate 
angles shows that the caps are close to elliptical. In Figure 5.6 we 
compare the function G(0) to the ratio H (0)/H(if/2), including points 
derived for intermediate angles. The differences between G(0) and our 
results are sufficiently small 10 percent) to be neglected. Thus we 

take over from Chapter 3 the theoretical distribution f(y) of pulsars 
in the variable y = (T VW) . f(y) peaks when W/T = 2l|(^>/2jt so 

= (^ 2 /9)(CM*/4Tt 2 ) 1/3 (i/R*) , (5.5.iOa) 


(2) = pi/3 (1) 

peak y peak 


(5-5-iob) 


for Criteria (1) and (2), respectively. We have evaluated these quanti- 
ties for a variety of neutron star models (Cohen and Cameron 1971 ; Baym, 
Pethic, and Sutherland I 97 I ; Leung and Wang I 97 I ; Clark et al . 1971 ) ; 
the results are presented in Figure 5 . 7 . 

The data in Table 5-1 have been taken from pulse profiles published 
by Manchester ( 1971 b) and Lyne, Smith, and Graham ( 1971 ) . We obtained the 
total width of the pulses by measuring the width at ten percent amplitude, 
and have included only those pulsars for which the noise at this level was 
small. In presenting the data, we employ "moving bins" (Turk 1972) to 
smooth out statistical fluctuations. In Figure 5-8, each data point repre- 
sents the number of pulsars in a bin centered on that point. The spacing 
between points is smaller than the bin size, so that each pulsar is counted 
more than once (seven times in 5 . 8 a and five times in 5 . 8 b). The error 
bars are just the square root of the number in each bin. The data peak at 
log (y) = 2.95 ± 0.10. From Figure 5-7 it is clear that all of the 
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neutron star models cover this region. The best fit occurs for Curves 
1 through the models of Cohen and Cameron (based on an earlier 
"stiffer" equation of state) are less satisfactory. To account for a 
distribution of neutron star masses we fold together theoretical curves 
corresponding to a range of neutron-star masses. The smooth curves in 
Figure 5*8 give the results for uniform distributions of log y of half- 
width A log(y) = 0, 0.10, and 0.20 around log(y ) = 2.95- The 
range of masses to which these correspond may be read off Figure 5.7, 
for the five distinct models of neutron-star structures. 


VI . Interpulses 

There are two tests of this model concerning those pulsars with 
interpulses. First, following the procedure outlined in Chapter 3 and 
the results of Appendix D, we may calculate the number of pulsars ex- 
pected to show interpulses by summing = P^/P^ over the 89 known 
pulsars (Terzian 1973) • Putting R^. = pR^, , 


P. 

x 


6.79 P' 1/2 *i /2 Ml l/6 T- X/ 3 ; 


( 5 . 6 . 1 ) 


the resulting sum of P^ depends on the neutron star model. Choosing 

those of Baym, Pethic and Sutherland (I 971 ), the number of pulsars 

expected to show interpulses lies between 1.6 (M = 1.4l M ) and 12.4 

(M = O.O925 M ) (taking p = 1) . This is illustrated as a function of 
* o 

the neutron star mass in Figure 5<9> and is in reasonable agreement with 
the existence of six pulsars with interpulses. 

The second test is a consequence of the fact that a pulsar showing 
two pulses must have © « 90 °, and thus the pulse width will satisfy 


W/T ^ 2¥j>/2rt. In other words, these pulsars must have y s y peak • 


65 



pulsars 0531+21, 0823+26, 095^08> and 1929*" 10, comparison of Figure 5.7 
and Table 5.1 shows that there are neutron star masses for which this 
condition is met. We have no reliable pulse width for PSR 0904+77 nor 
PSR 1055-21, which have periods 1.579 s and O.I97IS, respectively. For 
neutron star masses between 0.0925 and 1.4l M^, our model predicts 
pulse widths less than (320 msec - 43 msec) and ( 8l msec - 11 msec), 
respectively. [Terzian (1973) lists upper limits of 80 msec and 30 msec 
for the fifty percent widths of PSR 0904+77 and PSR 1055“21 ; these 
translate into widths at ten percent of peak of about 140 msec and 5 ^ msec, 
respectively, in agreement with our estimates.] 

VII . Pi scussion 

Although we have used the idealized model of dipole field lines and 
single particle motion to study the collection of non-relativistic gas 
in pulsar magnetospheres, it appears that the results are substantially 
independent of these assumptions. Material is trapped stably when the 
magnetic field line along which it moves cuts the surfaces of constant 
potential energy in such a way that there exists a point at lower poten- 
tial than neighboring points in either direction. As distance from the 
neutron star increases beyond R^, the equipotentials approach cylinders 
aligned with the rotation axis. Any closed field line which extends to 
this distance will have such points, regardless of its detailed shape. 

It appears that, for almost any field configuration, the field will trap 
gas starting somewhere near r = R^, and the field will be opened 
shortly beyond that radius. The magnetic torque on the neutron star 
will thus be that calculated in Chapter 4, with a braking index of 

n = 7/3. 


66 



The assumption of single-particle motion is good provided the term 
(p/m) W in the equation of motion dominates the term V||P, where p 
is the pressure. If d = x R^, is the characteristic size of a collec- 
tion region and |w| = y F , this requires 

03 2 R FB 2 xy » kT|,/m . (5-7-1) 

For parameters of the Crab pulsar, Eq . ( 5 . 7.1) implies that the model 
is acceptable if 

T|| « 10 10 ' 9 (m/m ) xy . (5*7.2) 

O q O 

Reasonable values for x and y are 0.1, so that we need T « 10°*^ K. 
We believe that this is not a very restrictive requirement, and infer 
that gas collected around stable points will not be forced out of those 
regions by its own pressure. 

We have also neglected the mirroring force 

| v ^ 2 if 1 V, t B . (5.7.3) 

This is much smaller than |(p/m) w| if 

kT^/rn « A ^ 2 xy , (5-7 -^) 

the same criterion satisfied by the parallel temperature in Eq. (5-7. 1 ). 

Containment of gas also requires that diffusion of plasma across 
field lines is slow. From Spitzer ( 1962 ), the velocity of transverse 
diffusion is approximately 

v ± D = 10 ' 2 ' 7 InA B ~ 2 T _1/2 V x (p/m) . (5. 7-5) 
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Taking the Coulomb logarithm InA to be 5, the time scale for transverse 
diffusion is about 


D 


T 

1 



10 


13.8 


2 

X KS 


10 


11.8 


sec 


(5-7.6) 


Apparently, transverse diffusion is not a problem. However, since the 
rate of mass loss from collection regions determines the post-glitch 
relaxation time (see Chapter k) , more detailed study is needed of these 
and related processes . 

From these discussions it seems that if there is a substantial amount 
of cool gas in pulsar magnetospheres, the mass and field configurations 
will be similar to those we have described. 
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TAELE CAPTION 


Table (5.1). The 25 pulsars included in our sample. The ten percent 
widths measured from the papers of Manchester (1971b) 
("RNM”) and Lyne, Smith, and Graham (1971) ("LS&G”) are 
listed under the radio frequencies (MHz) used for each 
observation. Widths which are underlined were judged 
the most reliable and those in parenthesis the least 
reliable. Our assessment of the "BEST" width takes 
this into consideration. The last two columns list 
the values of [W( sec) ]/[T( sec) ] , and log^y, y = (T'^'Vw)^. 
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Table 5.1 


PULSAR 

T( sec) 


TEN PERCENT 

PULSE 

WIDTH I 

'msec) 


W/T 

log y 



RNM 

— 

L S 

& G— 

— 






4i0 

1665 

151 

240 

4o8 

610 

"best" 



0329+54 

0.7145 

54 

5g 


61 

55 

53 

55 

.077 

2.32 

0450-17 

O.5498 





26 


26 

.047 

2.82 

0525+21 

3.7455 






225 

225 

.060 

2.0 6 

0531+21 

0.0331 

(4.7) 

1 





1 

.030 

4.03 

0628-28 

1 .2444 






145 

145 

.117 

1.80 

0809+74 

1 .2922 

86 

(io4) 

110 

95 



90 

.070 

2.24 

0818-13 

1 .2381 





3§ 


36 

o\ 

OJ 

O 

3.01 

0823+26 

0.5307 


12 


17 

15 

16 

16 

.030 

3.22 

0833-45 

0.0892 


3-5 





3-5 

•039 

3.51 

0834+06 

1.2738 

32 

(35) 

33 


33 


33 

.026 

3.10 

O95O+O8 

O.2531 

20 

20 

35 

19 

(31) 

21 

20 

• 079 

2.60 

1133+16 

1.1879 

4o 

33 

48 

43 


38 

4o 

.034 

2.90 

1237+25 

1 . 3824 

62 

54 

(68) 

62 

57 


60 

.043 

2.63 

1508+55 

0.7397 

2 7 

(29) 



26 


26 

•035 

3.00 

1604-03 

0.4218 





12 


19 

.045 

2.94 

1642-03 

0 . 3877 

6.5 

(11) 



8 


7 

.018 

3.76 

1706-16 

O.6531 

22 




(24) 


22 

.034 

3-07 

1818-04 

O.5981 

22 

(23) 





22 

.037 

3-02 

1911-04 

0 . 8259 

14 

16 



14 


i4 

.017 

3. 6o 

1919+21 

1 • 3373 

45 

(48) 


4i 



43 

.032 

2.92 

1929+10 

0.2265 

11 

11 



18 


13 

• 057 

2.91 

1933+16 

0 . 3587 

13 

8 





12 

•033 

3.25 

2016+28 

0 . 5580 

22 

22 



21 


22 

.039 

2.98 

2021+51 

0.5292 

29 

22 





25 

.047 

2.83 

2045-16 

1.9616 

93 

82 


(93) 

96 


90 

.046 

2.48 


TO 



FIGURE CAPTIONS 


Figure 5.1. Structure of the magnetosphere in the aligned case. Entire 
shaded section — collection region with negative energy. 
Heavy shading — collection region with inward radial force, 
(a) Emission cone if Criterion 1 is adopted, (b) Emission 
cone for Criterion 2. (c) Emission cone if field lines 

are dipolar out to the light cylinder. (We have used Crab 
pulsar parameters in this figure.) 


Figure 5.2. Left , the magnetosphere model of Goldreich and Julian (1969). 

Right , this model. The magnetic field configuration is sche- 
matic, but the scale (Crab pulsar parameters) is approximately 
correct, R^: R^: R^ « 100:10:1. 

Figure 5.3* Division of (T|, (3) -space according to (a) the number of zeros 
of Q(x) i n interval 0 < x < 1» ( b ) sign of e 
(clear area, negative everywhere; hatched area, positive 
around x = 1* becoming negative somewhere in the interval 
0 < x < 1; shaded area , negative around x = i> becoming 
positive and then negative again somewhere in the interval 
0 < x < i) : ( c ) the sign of f (same notation as Figure 3b) . 

Figure 5 .k. The roots x of Q,R, and S, plotted versus (3, for 

several values of Tj. Entire shaded region , positive radial 
force. Heavy shading , positive energy. The heavy curves 
are the stable points of single-particle motion. The 
features described in the text are illustrated here. 
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Figure 5.5. Polar caps, plotted as r = k(<|>) versus (j). Broken lines , 
aligned case. Solid lines, orthogonal case. In both cases 
the larger polar cap corresponds to Criterion 1 , the smaller 
one to Criterion 2 . 


Figure 5.6. Comparison of the function G(@) [smooth curve] with the 
ratio of minor to major axes derived in this model [ open 
circles, Criterion 1 ; filled circles, Criterion 2 ]. 


Figure 5.7. 


The quantities log (y^ 1 ^ ) and log (y^ ) for several 

l(j p6aK iU pcaK 

neutron star models: Curve 1 , Leung and Wang (1971) [their 
curve 1 ], Curve 2 , Leung and Wang [their curve 2 ]. Curve 
3 , Baym, Pethic, and Sutherland ( 1971 )- Curve 4 , Clark et al, 
( 1971 ) [their table 5]. Curve 5 , Cohen and Cameron (I971) . 


Figure 5 * 8 - Comparison of the observed distribution of pulsars in the 

2/3 2 

variable y = (T J /W) with the results of this model. 
See the text for an explanation of the symbols. 


Figure 5.9. The number of pulsars (out of the presently known 89 pulsars) 
expected to show interpulses as a function of the neutron 
star mass, for the neutron star models of Baym, Pethic, and 
Sutherland (1971). 


72 








74 




75 


Figure 5.3 
























Figure 5.5 
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Figure 5*6 
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Chapter 6 


DISCUSSION— CONCLUSIONS AND PROSPECTUS 


In the Goldreich-Julian (19 69) steady-state model, charged particles 
in the closed magnetosphere are tied to magnetic field lines which coro- 
tate with the underlying neutron star. This follows from the equation of 
motion under the assumption that gravity and inertia are negligible com- 
pared with the electromagnetic forces. In our model on the other hand, 
gravitational and centrifugal forces are assumed to dominate in the region 
of closed field lines, but the plasma still is tied to corotating field 
lines. Thus it is reasonable to inquire as to the consistency of our 
calculation of microscopic stability. 

The electric field which drives corotation is derived from a charge 
density which may be found from Poisson's equation: 


ejn -n_| ~ E/4nr ; (6.1) 

r is a characteristic length in the magnetosphere. In a magnetic field 

■4 — * 

this produces an E x B drift 


V 

EXB 


cE/B~4jtcer n -n_ 



( 6 . 2 ) 


so that the condition of corotation, v ~ cor, implies that the charge 

exb 

density is 


n + -n ~ cnB/4jtec . 


(6.3) 


This is of course the same as that given by Goldreich and Julian [Eq. 

( 2 . 12 )]. 

The gravitational force per unit volume is 
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(6.4) 


f ~ GM V m (n +n )r 2 , 
g * P v + 


the centrifugal force per unit volume is 


f ~ ( n +n ) m cd r , 
c v p ’ 


(6-5) 


and the electromagnetic force per unit volume is 


f ~ n -n evB/c . 
em + 


( 6 . 6 ) 


Thus gravity dominates electromagnetic fo-rces if 


n +n 

+ - 


n -n 
+ - 


» 


/ r WeB ) 

Kb j VW 


(6.7) 


and centrifugal force dominates electromagnetic forces if 


n + +n 



» 


eB 


m cjdc 

P 


(6.8) 


[The criterion (6.8) is the same as that given by Ruderman (1972) in order 
that the Alfv^n cylinder be substantially inside the light cylinder.] Sub- 
stituting the charge density required for corotation from Eq . ( 6 . 3 ); the 
mass density must satisfy 


m p (n + +n_) » B 2 /4jtc 2 . (6.9) 

Evaluated at the force-balance radius for the Crab pulsar, the right 
hand side of Eq . ( 6 . 9 ) is about 10 g cm , A mass of lO^ 1 '^ g dumped 

from Rp B is appropriate for a glitch in the Crab; if contained in a 

2 -3 

volume V = n(Rj, B /5) 2n this has a density about 1 g cm . Thus 

it appears that the kind of mass densities involved in magnetospheric 
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activity insure the consistency of our assumptions about microscopic 
stability . 

The criteria which we have used to decide between macroscopic sta- 
bility and instability were based on the assumption that so much mass 
will collect that the magnetic field cannot support it against gravity 
or contain it against centrifugal force. We have seen that microscopic 
stability can occur only for radii comparable with or greater than the 
force-balance distance, so we will concern ourselves with the effects of 
centrifugal force . The magnetic force which contains plasma is 

f m = V(B 2 /8jt) ~ B 2 /8«r . (6.10) 


We may define a critical mass density p c by equating and f m , 


p c = (B 2 /8«c 2 )(R L /r) 2 ; 


( 6 . 11 ) 


densities which exceed this cannot be held in corotation except in the 
force-balance region. 

The maximum plasma moment of inertia which can be contained by the 
magnetic field outside of the region of cancellation of gravitational and 
centrifugal forces is 

c \ 

'max = I p/ “tar 2 dr - B^/6a>% B 3 • (6.12) 

*i'B 
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Evaluated for Crab pulsar parameters, I ~ 10 g cm , so that 

max 

—12 

I /i ^ 10 . This is three or four orders of magnitude too small to 

max * 

account for the glitches, but is suggestive of the noise seen in the 
Crab timing data. We infer that non-relativistic gas at densities greatly 
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exceeding p c collects near R^, , and that when densities approaching 
build up outside they result in macroscopic instability. The 

resulting changes in the magnetic field configuration may be related to 
the restless behavior of the Crab pulsar. 

Lehnert (i 960 ), in a paper on the confinement of rotating plasmas 
by a magnetic field, notes that "in a theory on the origin of the solar 
system, Alfv^n (1954) considers a cosmic cloud of ionized gas, rotating 
around a central mass M c in a magnetic dipole field which has a dipole 
moment coinciding with the axis of rotation. In the theory is introduced 
a dividing surface. At the inner side of this surface, the parts of the 
cloud are attracted and fall down to the central body. Outside of the 
same surface, the centrifugal force brings matter out along the magnetic 
field lines down to the equitorial plane where planets or satellites may 
be formed. Alfv&i determines the dividing surface from the balance 
between the gravitational and centrifugal forces ." His expression for 
the dividing surface, when rewritten in spherical polar coordinates, is 

r 3 sin 2 e = 2 G U q /3cd 2 ; ( 6 . 13 ) 

using Eqs . (5.2.4), (5.2.7)> (5-2.ll), and (5.2.25), this becomes 

X 4 = 2/3 p 3 - (6.14) 

This is identical with our Eq. (5.3-4), which defines the boundary of the 
microscopically stable region in the aligned case. 

Lehnert further states "In connection with this problem should also 
be pointed out that a weak magnetic field ... is expanded radially by 
the centrifugal force. The present deductions merely serve as an illus- 
tration to the situation where the magnetic field is strong enough to 
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balance this force." This is analogous to our assertion that outside of 

the force-balance region centrifugal force •will convert closed field 

lines into open field lines, and that the resulting magnetic field will 

-2 

fall off as B«r 


The pulsar model at which we have arrived combines the magnetospheric 
structure laid out in the preceding chapters with the radiation mecha- 
nisms described by Sturrock (197O, 1971a). The modifications of Sturrock's 
work which arise in the change from R^ = to R^ = R^ have been worked 
out in detail by Turk ( 1972 ) . A discussion of the results has been given 
by Roberts, Sturrock, and Turk (1973), and we will briefly mention some 
of their conclusions. 

Table 6.1 sums up some of the comparisons between a model where the 
radiation arises at the light cylinder and the two polar cap models — 
first, that for which R = R^, and second that for which R Y = R^. 
Concerning the braking index, the first two give n = 3 , which is wrong. 

The revised model gives n = 7/3 = 2.33, which is consistent with the 
data. The fact that some pulsars have an interpulse which is off-center 
is hard to understand in terms of the light cylinder model (this is 
because multipole moments higher than dipole will have fallen off faster 
at the light cylinder, unless the contribution to the magnetic field of 
currents flowing in the open magnetosphere is significant; Sturrock 1971 a) . 
However, it follows naturally from either polar cap model. Further, the 
revised polar cap model enables one to understand the number of pulsars 


which show an interpulse. 

Timing noise is not a natural consequence of either of the first two 
models, but it is a readily understandable property of the third. Since 
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radiation is produced near the surface and the magnetospheric configuration 


is determined near the force-balance radius, it is easy to understand the 
fact that the pulse-shape of the Crab pulsar did not change around the time 
of the September 1969 glitch. On the other hand, we do not think it is 
easy to understand a large glitch, such as have been observed in the Vela 
timing data, as being a magnetospheric phenomenon. Concerning the period- 
pulse-width relation, there is no clear indication of what this should be 
for the first model; we know that the second model gives one which is 
incorrect, but it appears that the third model gives a relation which is 
compatible with observation. There is growing evidence, however, that 
there may be two distinct classes of pulsars, and the W-T distribution 
should be re-examined in this light. 

None of these models offers an obvious explanation of the drifting 
subpulse phenomenon. However, we have mentioned one possibility in 
Chapter 5 , and the presence of non-relativistic gas in the magnetosphere 
suggests another— rthat Alfv^n-wave oscillations might be involved. These 
ideas should be explored more fully. 

The radiation properties of the three models are also summarized in 
Table 6.1. Little is known about the electrodynamics of the light cylinder 
model, so it is hard to make comparisons with the data. Looking at the 
two polar cap models, we see that either explains the period-age distri- 
bution, and that either can account for the RF properties. According to 
Turk, Sturrock, and Petrosian (1973)? the force-balance model can better 
explain the optical properties. Either can explain the x-ray properties 
and the particle flux into the Crab nebula. 
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TABLE 6.1 


MAGNETOSPHERIC STRUCTURE 
AND RADIATION MECHANISMS 


LIGHT POLAR CAP POLAR CAP 
CYLINDER \ = \ \ = r fb 

1. BRAKING INDEX X X V 

2. INTERPULSE, OFF-CENTER X V V 

3. TIMING NOISE, CRAB GLITCH X X y/ 

k. GLITCH (BIG, VELA-TYPE) ? ? ? 

5. PERIOD-PULSE-WIDTH DISTRIBUTION ? X ^ 

6. DRIFTING SUB-PULSE STRUCTURE XX? 

7. PERIOD-AGE DISTRIBUTION ? y/ V 

8. RF PROPERTIES ? V y/ 

9. OPTICAL PROPERTIES ? X y/ 

10. X-RAY PROPERTIES ? y/ V 

11 . PARTICLE FLUX INTO CRAB NEBULA ? y/ y/ 

Author's evaluation of properties of models: \/ indicates apparent agree- 

ment with observational data; X indicates apparent disagreement; ? 
indicates uncertainty. 
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Appendix A 


POLAR CAPS AND EMISSION CONES 


In a coordinate system where the magnetic axis is in the z-direction 
and the rotation axis lies in the x-z plane at polar angle 0, dipole 
field lines are the surfaces 


-1 2 „ -1 
r sin 0 = r = constant , 
o 


(A.l) 


and a cylinder of radius R^ aligned with cn is the surface 


r s 


iny = Ry 


(A. 2) 


where y is the angle between r and to, 


cosy = cos0 cos© + sinQ sin® cos(j) . 


(A. 3) 


For a given (<j);®), these two surfaces intersect at ( r, 0) determined by 


2 

r Q sin 0 siny = R^. 


(A. 4) 


The first field line to intersect the cylinder (the one which is tangent 

to it) does so at an angle 0 such that r is a minimum. The condi- 

c o 

tion 3r /90 =0 is 
o 

sin0[2cos0 siny + sin0( cosy/siny) ( sin0 cos© - cos0 sin© cos(j))] = 0 . 

(A. 5) 

(The solution 0=0 corresponds to a maximum of r Q • ) 

Tracing these grazing field lines back to the star's surface r = R^, 
they form two polar caps, at colatitudes 0^ and it - 0 f , where 

—1 P “1 2 —1 2 

R„ sin 0 = r sin 0 = R sin 0 siny( 0 ) , (A . 6 ) 

* rc cY c 'c' 
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and two emission cones consisting of the tangents to the field lines 
where they cut the surface. The half-angle \Ji of a cone is given by 

tani|> = (3/2)sin0 r {[l - sin0 r ]/[l - ( 3/2) sin 2 ^] } , (A. 7) 

which is well approximated (to order [W/T]) by 

t = (3/2) sin0 r . (A. 8) 

The solution of Equation (A. 5) for general ((();©) is involved, and 
we shall content ourselves with the major and minor axes of the emission 
cone for arbitrary ©, and the complete solution for © = rt/2. 

y-z Plane ((}) = n/2) , Major Axis ; 

With cos(j> =0, Eq . (A. 5) is 

cos0(2 sin^® + 3 sin^0 cos^0) = 0 , (A. 9 ) 

the only real solution of which is 0 c = jt/2, siriy = 1, so the half- 
angle of the emission cone in the direction of increasing longitude 
(referred to the rotation axis) is 

= (3/2) (VV 1/2 * (A ' 10) 


x-z Plane (o = 0) , Minor Axis : 

With cos(j) =1, Eq . (A. 3) implies 

Y = 9 - © , (A .11) 

so Eq. (A. 5 ) reads 
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(A. 12) 


2 

2 sinQ cos0 sin( 0-©) + sin 0 cos(0-®) = 0 . 

From this we find the relation 

tan® = (3 tan0 c )/(2 - tan 2 0 c ) . (A.13) 

(The expression given in Paper I contains a typographical error.) Noting 
that 0 c S rt/2, we invert Eq . (A. 13 ) to find 

0^ = arctan{[ -3~(9 4 ' 8 tan 2 ®) 1//2 ]/[ 2 tan®]} . (A.l4) 

Thus the half-angle in the direction of increasing latitude (referred 
to the rotation axis) is 

♦ b (0) = (3/2)(VV 1/2 G(@) ' (A ’ 15) 

where G is defined by 

\/0 

G( ®) = sin© sin ( 0 -®) , (A.l6) 

c o 

with 0 c determined by Eq. (A.l4). 

Polar Caps in the Orthogonal Case: 

When ®= jt/2, Eq . (A. 5 ) has 

P P 

sin0 cos0(2 - 3 sin^0 cos (J>) = 0 . (A. 17 ) 

2 

One solution is sin 0=1, the other 

c ’ 

sin 2 0 c = 2/(3 cos 2 lj)) , (A. 18) 

p 

for which we must have cos^fj) ^ 2/3* We ma y verify that Eq . (A.l8) indeed 
corresponds to the minimum r for (1) in this range. Defining 

O ! 
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(A. 19) 


F(cos 2 <j>) =[r o (sin 2 e c = l)J - |r o ( sin 2 0 c = 2/3 cos 2 (j))J , 

we find F(2/3) = 0, F '(x)=( 1/2) ( 1-x) “ 3/2 -(l/2)3 3/2 . Thus F'(x s 2/3) 

S 0, so that F(x S 2/3) ^ 0, QED. The resulting polar cap has the shape 


( r */ r y ) 1/2 2 1/2 3 3/ Vcos(j) , cos 2 (f) ;> 2/3 , 


sin0 ss 
r 


(A. 20) 


.(R*/^) 172 sin 1/2 (|) 


Q 

, cos (|) ^ 2/3 • 


The simplified model of Chapter 4 has = R^, and we ma y com- 
pare the resulting polar cap to those found in the more complete model 
of Chapter 5 * I n Figure A.l we have shown the scaled polar cap size 
x(<|>) predicted by Eq . (A. 20), along with those from Chapter 5> for the 
orthogonal case. The qualitative similarity occurs because the equipo- 
tentials of Chapter 5 are nearly cylindrical in shape, and the quanti- 
tative differences because they are perturbed by the gravitational poten- 
tial . 


Case of General (j) and 0 ; 

For completeness, we note that for arbitrary ( <J) ; ©) we may transform 
Eq . (A. 5 ) into a cubic equation in x = tan0: 


where 


3 2 

a x J + b x + cx+d = 0, 
a = sin® cos ® cos(f) , 

Q p 

b = 3 - sin^®(l + 3 cos <j)) , 

c = -5 a , 

2 

d = 2 sin © . 


(A. 21) 


(A. 22) 
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FIGURE CAPTION 


Figure A.l. Polar cap in the orthogonal case, plotted as r = hC*! 5 ) 

versus <J>. Solid lines , the model of Chapter 4, where 
= Rpg- Broken lines , the model of Chapter 5 . The 
larger cap corresponds to Criterion 1, the smaller one 
to Criterion 2. 
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Appendix B 


THE PERIOD-PULSE -WIDTH DISTRIBUTION 


Here we derive the period-pulse-width distribution expected from 
two emission cones located at opposite magnetic poles. Measuring angles 
from the rotation axis, let the cones be located at colatitudes © and 
« - 0, and have longitude and latitude half-angles ^ and respec- 

tively. Assuming the shape of the cones to be approximately elliptical, 
and that i) and are sufficiently small that plane geometry may 

be used, the pulse width W is related to the rotation period T by 

T 

2 it sin 


_ 

♦b 2 


2-.1/2 


(B.l) 


and a similar expression with © -• it - ® for the other pole. Here (|) 
is the colatitude of the line of sight of the observer. If neither 
expression is found to be real, radiation is not detectable by the 
observer. This geometry is illustrated in Figure B.l. 

The probability that W is in the range W to W + dW, given by 
twice the probability for one cone, is P w (W,T)dW, where 


P W (W,T) = 2 j* \ sin© d© f | sin(j) d^> 8 |w - [l - 


> 2 - 11/2 


[real integrand] 

We have assumed that both the locations © of the emission cones and 
the angle of the line of sight (j) are randomly distributed on (0,n). 
To transform the delta function we use the well-known rule 


(B.2) 


«[*(*)] = ^ 


df 

dx 


n 


1-1 


lx=x 


5(x-x n ) 


(B.3) 


95 



where 


f(x ) = 0 . 
v n ' 


(B.k) 


Defining 


°( 0 ’^ = * - ^n5f l 1 - ^2~\ /2 


(B.5) 


we have 


so T \ r 

5<f> it sin<J) 


! .ifca 
K 2 


21 


■ 1/2 




- 2 . 


T ncostj) r 


2 

it sinnj) L i|f 


i - 


2-1 


1/2 


(B.6) 


Here <J) are the solution of G = 0. i.e., of 
n 


W = 


h 


it sind) 1 1 2 

T n L ib. 


< 4 v®> 


2 1/2 


(B.T) 


Since we assume tj and ijr^ are small, we may put sinlj)^ ®* sin®, and 


thus 


^n = ® ± *b \ 1 " 


w 2 * 2 - 2 


1 


sin ® 


T V 


(b.8) 


and (B.6) becomes 


3G 


T 2 *„ 2 |l - 


M n 


= ± 


T V 


lb, W it 2 sin 2 ® 
T b 


3 

+ W sin ©cos® 


(b.9) 


To order (W/t) the first term dominates, so both solutions to G = 0 
are the same, and we may easily do the delta-function integration over 



P W (W,T) = T ‘ 2« J 


^ A sin^®(i)r b /^)d© 


2 1 

' - sin® 


W 2 / 


( B.10) 


where 


'rt 

i 2 


®k = 


sin 


-1 



*W 


atfW / ' irW 


s l 


< l 


( B. 11) 


We have used the invariance of the integrand under © -♦ jt - ®. 
In our model 

* _ 3 R 1 / 2 r - 1/2 
H - 2 R * \ 


♦ b = ty, G (®) 


(B.12) 


We 


may parametrize R^T) b y R 0 and 11 by defining 


»V - V 11 ; 


(b.13) 


if we further define 


9 R * 


k = 


4jt 2 R 

o 


(B.U) 


y = T 2_T1 W -2 , 


(B.15) 


so that 




V 


2 2 

Wjt 


= ky , 


(B.16) 
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then 


P W (W,T) = 2*t “ 1 j 


sin^© G(©)d© 
_ ^ 


o vky - sin 2 © 


It is convenient to transform to the probability that 3 
range y to y+dy. Using the identity 


|P w (W,T)dW| = |P y (y,T)dy| 


we find the central result 


P (y,T) = i T~^ 2 y _3/2 L(ky) , 


where we have defined 


L(x) = 2 « j* 
o 


sin^® G(©)d© 


i. 


x - sin ® 


with 


Si- 


ft 

2 


, X 2: 1 


-1 


sin x , x < 1 


Calculation of the function L(x) and its integrals is 
in detail in Appendix C. 


(B.17) 

is in the 


(B.l 8 ) 


(B.19) 


( B.20) 


(B.21) 


discussed 
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FIGURE CAPTION 


Figure B.l. The geometry of the polar cap model, showing only one 
magnetic pole. 
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Appendix C 


THE FUNCTION L(x) 


The function L(x) is defined by 


L(x) = 2 I t| m Sln4 . Q - G ( @ ) dQ 


where 


I — ' 

yx-sin © 


, x s 1 


0 = 
m 


-1 


sin > x < 1 


(C.l) 


(C.2) 


We first demonstrate that L is integrable near the singularity at 
x = 1. For x s 1 we use z = sin® to rewrite L as 


.1 . -1 

_ / \ „ C z G(sxn z)dz , sr. 

L(x) = 2jt I t= ~ (x ^ 1) , (C.3) 

o y(l-z 2 ) (x-z 2 ) 
and for x ^ 1 we employ sin® = z{x to write 

L(x) . 2* f 1 (is 1) . (C.4) 

o \/( l/x-z 2 )(l-z 2 )' 

The complete elliptic integral of the first kind is defined (Abramowitz 
and Stegun 1970 , 17-3.1) as 


K ( m ) = 


r 


[(l-t 2 )(l-mt 2 )]” 1/2 dt 


(C-5) 
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so 


J [(l-t 2 )(m-t 2 )]“ 1/2 dt = K (i) . 


As z and G are less than one, from (C.3) and (C.4) we find 


L(x) <; 2n x" 1/2 (x s l) , 


and 


L(x) 2 it x^ K(x) (x ^ 1) . 

The asymptotic form of K(m) is (Abramowitz and Stegun 17.3*26) 


K(m) 


m -♦ 1 


. l , 16 
t 2 ln l^m > 


so putting 


we have 


x = i ± e (x £ i) 


l(x) ^ 2K(i+e)' 1/2 | in(i6/e) , (x = i+e > 1 ) , 


L(x) <; 2n (1-e) 2 ~ m(l6/£) , (x = 1-t < 1) . 


It is evident that we need consider only one case, so we examine 


— 1 /? 

I L(x)dx ^ jt I (1+e) ln(l6/g.)dg . 


(c.6) 

(c.7) 

(c.8) 

( c *9) 

(c.io) 

(C.ll) 

(C.12) 

(c.13) 
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For small £ this is dominated by -jtg. In £, which shows that L(x) 
is integrable across x = 1. 

In order to tabulate the function L and its integrals we use the 
forms (C.3) and (C.4) . Subtracting off the singularities at z = 1, 


and 


L(x) = 2n 


g^ 2 ) = 


j* 1 g > ( z)dz + 2rt 
O 

4„, . -i v 

z G( sin z) 

\l( l-z 2 )(x-z 2 )' 


G( Tt/2) 

\/ 2 (l-z)(x-l)' 


(x ;> 1) (C.l4) 


L(x) = 2 w x 3/2 j g< (z)dz + 2rt x 3/2 G(sin , 

o 

(X g i) (c.15) 

, . z\j[sin 1 (z\/x)] Gfsin 1 \/x) 
g,(z) = . S — ■ 7 — • 

\/(1-z 2 )(1A-z 2 ) 

The remaining numerical integrals may be handled easily by standard 
methods. The resulting function L(x) may be integrated numerically, 
as we have shown that the contribution from the singularity at x = 1 
may be made arbitrarily small. The precise integral required is [see 
Eq. (3.16)] 

f z 2 

A( z 1 ,z 2 ) = J z” 3/2 L(z)dz . (C.l6) 

Z 1 

From Eqs . (C.3) and (C.4) it is easily demonstrated that 


t / \ „ 2 -1/2 

L(x) ^ n x 

X » 1 


(c.17) 
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and 


L( x ) 


S3 


2 

jt 


2 


x 


X « 1 


(c. 18) 


so that defined by 



(c.19) 


is finite (X = 7*772). The function A is recorded in Table C for steps 
of 0.1 in log^^z), and is illustrated in Figure C.l. 
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TABLE CAPTION 


Table C. The function ^{z^,z^) , as defined by Eq. (C.l6), for steps 
of 0.1 in log (z). The value of A for each interval is 
the number in column 3 times 10 to the power entered in 
column 4 . 
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TABLE C.l 


log z ± 

log z 2 


-2.0 

-1.9 

.131302 

-1.9 

-1.8 

.185621 

-1.8 

-1.7 

.262470 

-1.7 

-1.6 

.371237 

-1 .6 

-1.5 

.525261 

-1.5 

-1 .4 

.743518 

-1 .4 

-1.3 

.105307 

-1.3 

-1.2 

.149258 

-1.2 

-1.1 

.211752 

-1 .1 

-1.0 

.300781 

-1 .0 

-0.9 

.427929 

-0.9 

-0.8 

.610122 

-0.8 

-0.7 

.872385 

-0.7 

-0 .6 

.125231 

-0.6 

-0.5 

.180777 

-0.5 

-0.4 

.263121 

-0.4 

-0.3 

.387964 

-0.3 

-0.2 

.58^078 

-0.2 

-0.1 

.925510 

-0.1 

0.0 

.182218 

0.0 

0.1 

.177488 

0.1 

0.2 

.945468 

0.2 

0.3 

.647181 

0-3 

0.4 

.470388 

0.4 

0.5 

. 352098 

0.5 

0.6 

.268172 

0.6 

0.7 

.206565 

0.7 

0.8 

.160350 

0.8 

0.9 

.125165 

0.9 

1.0 

.980978 

1.0 

1 .1 
1.2 

1.1 

1.2 
1-3 

.77H77 

.607641 

.479629 

1-3 

1.4 

. .379101 

1.4 

1.5 

.299959 

1.5 

1 .6 

.237536 

1 .6 

1.7 

.188224 

1.7 

1.8 

.149225 

1.8 

1-9 

.11835^ 

1-9 

2.0 

.938993 

2.0 

2.1 

.745159 

2.1 

2.2 

.591454 

2.2 

2.3 

.469527 

2.3 

2.4 

.372782 

2.4 

2.5 

.295999 

2-5 

2.6 

.235050 

2.6 

2.7 

.186663 

2.7 

2.8 

.148244 



SUMS OF BINS 


. 372073 


.383792 



.098719 


.036114 


.013962 


.005497 


002179 


.000866 
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figure caption 


Figure C. 


Histogram of the functioe 4 ( 3 ^, as defined by 
Eq. (C.16). 
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Appendix D 
INTERPULSES 


We first consider the fraction of pulsars which are observed. If 
the (+) magnetic pole is at colatitude © and the (-) pole at colatitude 
« - ®, the observer detects radiation from the (+) pole if the line of 
sight <j) satisfies 


.® - * b < 4 < ® + t b > C 0 - 1 ) 

and he sees the (-) pole if 

lt-M b < <j> < *-0% . (D.2) 


These regions are shown in Figure D.l. Assuming random orientations of 
© and <J), the probabilities that he sees the (+) or (-) poles are 


®HL 

f it 1 f b 1 

— sin® d® l — sin(j) d<j> , 

®-t b 


(D.3) 


and 


jt-®# 

J rt 1 C b i 

— sin® d® I — sin(j) d(j) . 


(D.4) 


Clearly P + = P , and we evaluate only P + : 


1 

P + = tj; i sin®[-cos(®ti|f b ) + cos(@~t b )]d® , 


f n/2 

■> + = J 


2 

sin ® sinto_ d® . 

D 


(D.5a) 

(D.5b) 
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1/2 


Now sini|r = (3/2)(R^/R ) G(@), so 


P = 


= (3/2)(R^/R y ) 1/2 j* sin 2 © G(0) d© 


and the total probability that we see at least one pulse is 


1/P p 

p ! = 2P + = 3(R*/R y ) l sin ® G( ©) d© 


By numerical integration 


it/ 2 


sin^® G( ®) d® = 0.620958 , 


so 


P X =1.8629 (R*/Ry) 


1/2 


[Of course, we could also obtain P y from 


P l(T) = p y (y,T) dy 


which becomes ( see Appendix B) , 


P 1 (T) = 1/2 t" 1 ^ 2 k 1/2 X 


Using the definitions of k and R , we find that 

o 


t -T 1/2 fc l/2 = (3/ 2j t)(R^/R y ) 1/2 , 


so that Eq. (D.ll) becomes 


(D.6) 


(D.T) 


(d.8) 


(D.9) 


( D.10) 


(D.ll) 


(D.12) 
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(D.13) 


p x (t) = (3X/4rt)(R*/R y ) 1/2 . 

The numerical value of \ = 7-772, so that 3X/^Jt = I. 855 . The differ- 
ence between Eqs. (D. 9 ) and (D. 13 ) is about 0 .4$, perfectly satisfac- 
tory agreement considering the approximations used in the derivation of 

Eq. (B.19)]. 

We now consider the number of pulsars expected to show two pulses. 
The region of © — <)> space for which Eqs. (D.l) and (D.2) hold simul- 
taneously is illustrated in Figure D.l. When weighted with 
1/2 sin 0 1/2 sinlj), this area is the probability of a pulsar showing an 
interpulse, 



Clearly the two contributions are equal; evaluating the first gives 


i.i f 

2 2J 


it /2 


sin(j )[2 sin(|) sinty^-^) - 2 cos({) cos\|i h (-^) ] 


b K 2' 


«/ 2-<r b (§) 


(D.15*) 


P 2 "I Vf sin Vl) '^b 2( l ) = ‘°- 4 33 (V R Y ) ■ (D-^b) 

[This agrees with the observation that the weighting factor is very 
nearly 1/4 in the region of integration, and that the area of the square 
is (2V^) 2 = 2* b 2 (f).] 
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FIGURE CAPTION 


Figure D.l. The regions of (®,<j)) -space. Entire hatched area , observer 
detects at least one pulse. Cross-hatched area, observer 
detects both pulses . 
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n 3 


Figure D.l 
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